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ABSTRACT 



A deterministic simulation using a model of human dynamic 
orientation was written to optimize the parameters of the 
motion base control system for a six-degree-of-freedom flight 
simulator. An experiment requiring pilots to rate different 
levels of motion fidelity during a basic flight task pro- 
vided a data base for validation of the simulation. Ratings 
between subjects for linear, rotational, and combined motion 
cues were inconsistent due, in part, to the subjects' lack 
of experience in the F-15 aircraft and proficiency in high 
performance aircraft. The coefficient of concordance am.ong 
subjects for the three ratings were .4483, .4835, and .5914, 
respectively. Comparison of simulation results with ex- 
perimental data yielded positive correlations as high as 
.5138. Response of the simulation to changing wash-out 
filter parameters was investigated and found to be adaptable 
to experimental optimization methods. 
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I . BACKGROUND 



A- GENERAL 

The modern, complex flight simulator has assumed an 
expanded role in aviation research and development as well 
as in military and civilian training programs. In the face 
of the current energy crisis, manpower and budgetary cut- 
backs, and reduced flight operations, the Department of 
Defense (D.O.D.) has recognized the benefits of judiciously 
substituting simulator "flight” time for aircraft hours. 
Supplementing existing training programs with instruction 
in flight crew trainers has long been an effective means of 
training safe, professional aviators. D.O.D. has establish- 
ed an objective of reducing flight operations 25% by the end 
of Fiscal Year (FY) 1981. The Chief of Naval Operations 
(CNO) has recently authorized the substitution of simulator 
time for ten percent of total annual flight time requirements 
(25 percent of the annual instrument flight time) . Reference 
1 states the CNO's policy: 

"As additional simulators become available and more 
is learned on the 'transfer of learning' gained through 
the use of simulators, this program will be expanded." 

Increased emphasis on the use of synthetic training 
devices has accelerated the procurement and funding of 
special part-task trainers such as the Navy's Air Combat 
Maneuvering device and Universal Night Carrier Landing 
Simulator. Total D.O.D. FY77 procurement funding for train- 
ing devices is 298.7 million dollars as compared to 88.5 
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million in FY74. Research and development funding during the 
same period has nearly doubled [Ref. 2]. The development and 
deployment of modern combat aircraft today is paralleled by 
the production of high-fidelity weapon system trainers. Ad- 
vanced weapon platforms such as the F-14, S-3A, P3-C, SH-2F, 
and the SH-3H have (or will) become operational concurrently 
with advanced, full task, motion base trainers. 

The following section reviews the current "state of the 
art" of the flight simulator as a tool in the training en- 
vironment, research and development community, and in the 
field of accident investigation. It references selected 
studies and events that demonstrate the devices* capabilities 
and effectiveness. 

1 . Training 

Understanding the benefits of flight crew trainers is 
essential to their cost effective use within comprehensive 
training programs. Efficient use of simulators reduces the 
expense of fuel, weapons, support equipment, and manhours. 

It also reduces the competition for hardware between training, 
operational, and maintenance departments. Alloted training 
hours can be devoted entirely to the task at hand, eliminating 
time for such activities as preflight, refueling, and clear- 
ance acceptance. As well as reducing accidents, it provides 
a safer environment for practice of hazardous or emergency 
flight operations. The advanced system incorporates fea- 
tures such as "freezing" the problem for real-time critique, 

<?> 

replaying the maneuver allowing the student to objectively 
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review it while fresh in his mind, and computation of purely 
objective performance measures using computer based 
algorithms. (See Ref. 3 and 4 for examples.) 

Instruction in training devices provides better control 
of psychological factors such as stress and workload levels. 
Use of a select, limited number of experienced instructors 
enhances standardization of undergraduate and "fleet" pilots. 

The Naval Air Training Command, in 1972, began evaluation 
of a computer generated imagery system which complemented the 
existing TA-4 2F-90 flight simulator. The device provided 
visual scenarios such as aircraft carrier landings, a bomb- 
ing range complete with scoring system, and carrier catapult 
launches. The resultant success in the Familiarization Stage 
of the advanced student flight training syllabus was expressed 
in the following statements [Ref. 5] : 

"The major area of benefit has been the practice 
provided in VFR procedures and techniques. The addition 
of a carrier visual presentation has made student in- 
struction in CQ [carrier qualification] possible with 
much the same results as those experienced in FAM stage. 

The real benefit of the simulator however has proven to 
be in the weapons delivery phase of the syllabus. Bomb 
instructors have noticed a significant improvement in 
the students' ability to fly the pattern and develop the 
principles of bombing as a result of incorporating simu- 
lator flights into the weapons stage. ...After a short 
period of instruction, it was evident that students who 
had had several hops in the simulator were, to a signi- 
ficant degree, outperforming those who had not." 

"The real key to the value of the simulation is the 
hit-spotting program which provides the instructor with 
real time readouts of release altitude, release mach 
number, dive angle, and hit position. ...The instructor 
may also introduce any of a multitude of factors affect- 
ing mission completion including adverse weather, air- 
craft malfunctions, variable winds and many others 
affecting flight safety. ...The system has enabled the 
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instructors to provide the student aviator with a high 
degree of realism and a method of evaluation which, 
until now, has been unavailable even in the aircraft." 

2 . Research and Development 

The use of cockpit simulators in the research, devel- 
opment, test, and evaluation phases of nev/ aircraft allows 
the aeronautical engineer to develop flight systems and pro- 
cedures without use of the prototype aircraft. Concurrent 
development of mockups provides for efficient aircrew station 
design and analysis from a system and human factors stand- 
point. Design of displays and controls is simplified with 
construction of full-scale prototype crew stations. Realis- 
tic, motion base simulators, driven by mathematical models 
derived from wind tunnel testing, allow test pilots to opti- 
mize aircraft performance capabilities and "handling qualities" 
specific to the mission. 

Development of the U.S. Air Force/Boeing Advanced Medium 
STOL (short takeoff and landing) Transport was aided 
extensively by the use of Boeing fixed-base and NASA motion 
base flight simulators. The following advanced aircraft 
design concepts were evaluated on the devices before 
incorporation into the prototypes [Ref. 6]. 

■ engine bleed air for leading edge flap boundary 
layer control (BLC) 

■ redundant mechanical and electrical flight controls 

■ aerial delivery procedures 

■ satisfactory handling qualities and criteria for 
engine-out STOL approaches 



13 



■ "conventional" piloting techniques for STOL speed 
and flight path control. 

An added benefit, recently used in the development of the 
F-14, is the capability to design, evaluate, and optimize 
specialized sxabsystems. The Automatic Carrier Landing Sys- 
tem (ACLS) of the F-14 "TOMCAT" was developed using a 
3-degree of freedom (pitch, roll, and heave) NASA cockpit 
simulator [Ref. 7] . After "proper validation" of the flight 
dynamics, optimization of the parameters of the ACLS ' s three 
main components began. Use of simulated flights in lieu of 
aircraft flights yielded "significant increases and flexibil- 
ity in the number of combinations [of parameters] that can 
be examined for aircraft subsystems which have many parameters 
to specify. " Other advantages of the program were soon 
realized : 

"1. Allowed seven iterations of SPN-42 [shipboard 
guidance equipment] parameter values to be examined in 
two to three days instead of the normal two to three 
weeks and several flights. 

2. Provided qualitative predictions of pilot com- 
ments on closed loop ACLS control characteristics and 
quantitative prediction of performance statistics, pilot 
AQR (ACLS Quality Rating) , and to a lesser extent, open 
and closed loop frequency responses. 

3. Saved $120,000 and $150,000 and approximately 
six months in the SPN-42 optimization phase of ACLS 
development. " 

"The impact of this new T&E capability is indicated 
by the inclusion of specific periods of simulator util- 
ization in the ACLS development plan for the S-3A airplane. 
Similar utilization will likely be specified for the F-18 
airplane. In addition, the same NASA-Ames Research Center 
R&D simulator used for the F-14A has been reprogrammed to 
attack long standing ACLS problems of the A-7E and A-6 
aircraft . " 
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3 . Accident Investigation 



Cockpit simulators have proven to be an indispensible 
tool for reconstruction of pre-accident histories. The pres- 
ent capability to use data from inflight recorders to derive 
flight profiles, environmental conditions, and signals with 
which to "drive" simulators is impressive. Current digital 
recorders can record up to 64 flight parameters per second. 
However considerable re-formatting of the data is then re- 
quired in order that simulators can use the data directly. 

A current NASA program is in pursuit of the capability to 
recreate pre-accident flight conditions using information 
from flight recorders, voice recorder data, radar tracks, and 
meteorological data. 

The investigation into the crash of an Eastern Air- 
lines B-727 in June, 1974, by the National Transportation 
Safety Board included evaluation of the wind shear and down- 
draft activity encountered by the aircraft [Ref. 8] . Four- 
teen experienced pilots flew a total of 54 approaches in a 
fixed-base B-727 cockpit simulator programmed with four 
different hypothesized wind models. Even with an a priori 
knowledge of the type and severity of the storm activity, 

18 flights "crashed", while only five flights resulted in 
placement of the aircraft near the runway threshold. Subse- 
quently constructed flight profiles (airspeed and altitude 
traces) closely resembled those from the B-727 's flight 
recorder, validating, to some extent, the wind shear models. 
At the time of the crash, the flight crew was thought to be 
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attempting visual contact with the runway environment. Seven 
of ten pilots who commented following the tests felt that 
their efforts to "go visual" at decision height delayed their 
recognition of the sink rate produced by the wind shear and 
downdraft conditions. 

A committee, chaired by the Director of Aviation 
Safety Programs, Naval Postgraduate School, was formed in 
1975 to investigate present use and potential of flight sim- 
ulators in aircraft accident prevention and investigation. 
Representatives from government, NASA, and the airline indus- 
try recommended further research into the use of existing 
inertial navigation equipment to record available flight 
accelerations in a form adaptable to simulators [Ref. 9] . 

The committee offered these recommendations; 

"Crash-protected digital flight recorders should 
be installed in high valued aircraft (with inertial 
navigators) . There should be a unified effort to 
develop and implement standardized data collection and 
processing procedures which would facilitate recreation 
of accident/incident conditions in flight simulators. 

Voice recorders, especially on multi-crew aircraft, 
should also be considered." 



B. NATURE OF THE PROBLEM 

1 . Pilot-Simulator Feedback Loop 

While accomplishing a certain task in a flight sim- 
ulator, the pilot responds to various visual, kinesthetic, 
audio, and vestibular feedback cues (see Fig. 1) . Feedback 
may be classified as either "intrinsic" (e.g., control "feel" 
or amplitude of movement necessary in the operation of a 
control) or "extrinsic" (the consequences or results of 
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Figure 1; Pilot-Simulator Feedback Loop 



control operation) . The problem encountered by a designer is 
one of cost-effectiveness, i.e., which cues are necessary and 
what degree of fidelity is required of each cue to ensure a 
positive transfer of training at the particular stage of 
learning for the given task. 

The requirement for motion fidelity in a flight sim- 
ulator has been an ongoing question for years. It cannot be 
debated that motion cues simulating airframe buffeting as 
well as in response to flight control inputs adds greatly to 
psychological fidelity, especially during early stages of 
training. Many research studies have compared fixed-base 
and motion base simulators measuring pilot performance, 
learning rates, and subjective pilot responses. VJhere 
motion cues play an important role in the pilot's feedback 
loop, the addition of motion cues creates a significant 
improvement in performance. For example, during instrument 
approaches, the inherent lag in cockpit instrumentation 
forces the pilot to rely on the vestibular sense to a higher 
degree. (A possible argximent for fixed-base instrument 
training, where the pilot is forced to depend on artificial 
cockpit references) . During visual carrier landings, the 
first sensation of rapid glide path deviation is typically 
through the vestibular and kinesthetic sensors, quickly 
confirmed by movement of the visual landing aid. Another 
scenario in which motion cues are important to the pilot is 
low-level, high speed terrain following flights. Where motion 
tends to degrade pilot performance, the addition of these cues 
allows a quicker transfer of learning to the aircraft. 
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2 . Physiology of Motion Detection 



Man perceives the sensation of motion through three 
main sensors. The eyes view the changing world around him. 

The kinesthetic sensors within his body respond to the rate 
of change of movement of his body and its members. The semi- 
circular canals and otolith organs are the main components 
of the highly specialized vestibular system of the inner 
ear. The inner ear canals react to changes in rotational 
velocity, i.e., angular acceleration. The utricles and 
saccules of the otolith system respond to linear accelerations. 

The present intent is not to analyze in depth the 
physiology of the vestibular system, but to review the mo- 
tion to which the organs respond. (A more comprehensive 
review of vestibular physiology and its function is available 
in Ref. 1 0.) 

The fluid-filled semicircular canals form a roughly 
orthogonal set of sensors in each inner ear. The endolymph 
fluid lags behind the movement of the canal walls when the 
head undergoes rotation. This relative motion tends to dis- 
place the cupula which obstructs an expanded portion of the 
canal called the ampulla. Deflection of the cupula excites 
sensory hair cells at its base causing an increase in the 
firing rate of the afferent nerve fibers. Orientation of 
the canals is such that those on the left side of the head 
are coplanar with those on the right side of the head, 
sensing motion about the same axes, but in the opposite 
direction, thereby acting as a single unit. It can be 
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assumed that the central nervous system responds to the 
difference of their respective afferent responses. 

The otolith system of the inner ear is comprised of 
two utricles and two saccules. Each is a mass of calciiam 
carbonate crystals suspended by sensory hair cells. Linear 
accelerations tend to displace the mass exciting the af- 
ferent nerve endings at the base of the hair cells. The 
utricle organs are oriented so as to be sensitive to the 
horizontal component of specific force. The saccule organs, 
approximately perpendicular to the utricles, sense vertical 
shear forces and are roughly one-half as sensitive. 

3 . Vestibular Modelling Research 

The need for development of a model of the human 
vestibular system can be illustrated best when one considers 
the advances made in modern aerospace technology. The 
high G environment of combat aircraft and space launch 
vehicles, the prolonged weightlessness of spacecrafts, and 
the motion experienced in large ships and tall buildings are 
products of twentieth century technology requiring the 
understanding of man's capacity to function in these 
environments . 

The "state of the art" of current modelling capabil- 
ities is best illustrated in a PhD thesis by Ormsby [Ref. 11]. 
At the MIT Man-Vehicle Laboratory, under a NASA grant. Dr. 
Ormsby developed a mathematical model based on the known 
physiology of the human vestibular apparatus. The signal in 
noise model is capable of predicting near threshold and 
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supra-threshold perceived motion and attitude in the pres- 
ence of random noise and spontaneous firing of the afferent 
sensors. It was then refined to be consistent with current 
neurophysical and psychophysical data. Included is a model 
of the central processor (the brain) and its interpretation 
of afferent dynamics with an a priori knowledge of the mo- 
tion stimulus. Since it was difficult to obtain data 
resulting only from vestibular cues, the model is said to 
also reflect responses resulting from tactile and kinesthetic 
motion cues. 

Separate models of the semicircular canals and the 
otolith organs are combined to yield the following estimates 
in the coordinate system of the head: 

■ perceived linear acceleration 

■ perceived rotational velocity 

■ perceived attitude (a unit vector in the perceived 
"down" direction) . 

The model has been successful in predicting illusions 
such as the dynamic elevator illusion in the absence of vi- 
sual inputs. It also has produced results in close agreement 
with data describing perceived lateral tilt angle as a func- 
tion of actual tilt angle in "IG" and "2G" environments. The 
pitch and roll sensations during a catapult launch were 
investigated and used to further refine the model. 

It must be emphasized that the version of the model, 
as used in this thesis, is strictly a vestibular model with- 
out visual, tactile, or kinesthetic inputs. The integration 
of vestibular and visual cues is the subject of an ongoing 
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research program at the MIT Man-Vehicle Laboratory 
[Refs. 12 and 13] . 

4 . Components of a Motion Base Control System 

The human vestibular system is sensitive, not to 
motion, but to the higher derivatives of motion such as 
linear and angular acceleration. Thus only the "rate of 
change" of motion need be simulated. The simulation prob- 
lem, then, is how to optimally simulate acceleration cues 
given fixed operating limits. (Hardware has inherent size, 
velocity, and inertial limitations, while computers have 
speed and accuracy limits.) A computer model of the air- 
craft produces acceleration vectors in response to control 
inputs and external factors such as wind and turbulence 
which are converted to simulator drive commands. Linear 
accelerations are simulated by tilting the cockpit cab 
(such as backward on application of thrust) at a rate be- 
low human rotational thresholds. Thus lateral and longitu- 
dinal acceleration cues can be simulated, but only the 
higher frequency vertical cues (such as turbulence and 
rapid pullups) can be simulated. After the cue has been 
transmitted, the cab is gradually returned tc its neutral 
position. This tendency, known as "wash-out", keeps the 
cab near its center position to maximize its operating 
envelope for subsequent cues. 

This section describes the main components of motion 
base control systems. The vast majority of motion base 
simulators have either "synergistic" platforms supported and 
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driven by linear hydraulic actuators or electro-mechanical 
drive systems with independent degrees of freedom. Digital 
computer programs for both system types are comprised of 
three main components : 

■ aircraft models which produce the particular 
flight dynamics and control feedback ("feel”) 
characteristics 

■ "wash-out" filters which remove low frequency 
accelerations which tend to generate large 
amplitude movements 

■ "drive" programs which transform the aircraft 
motion to simulator motion commands. 

This sequential process is shown in Figure 1. The motion 
of the aircraft about its six degrees of freedom is used to 
drive cockpit instrumentation, update navigational equipment, 
and drive external visual displays. The drive program 
consists of algorithms which compute: 

■ residual tilt which uses gravitational forces 
to simulate linear accelerations 

■ lead compensation for "quickening" the position 
commands to compensate for inertial lag in the 
machine ' s dynamics 

■ cab-to-drive mechanism transformations 

■ "soft" position, velocity, and acceleration 
limiting networks. 

Although many additional safety features and utility func- 
tions are present in most programs, these are the primary 



23 



motion generating components. (See Refs. 14 and 15 for more 
in-depth, theoretical descriptions of control systems.) 

5 . Simulator Fidelity 

Studies which evaluate motion control systems employ 
sophisticated means of calculating performance measures and 
compiling subjective pilot responses. Numerous studies have 
evaluated the handling quality differences between fixed 
and motion base simulators in an effort to understand the 
need for and the contribution of motion cues to the pilot. 
(See Refs. 15 and 16 for examples.) Improving stability 
characteristics for new aircraft or simulators often in- 
volves the analysis of pilot performance while system param- 
eters are systematically changed. Verbal responses are 
analyzed to discern which aspects of the motion cues are 
deficient. Design of an aircraft control system typically 
entails a series of compromises and tradeoffs yielding 
handling qualities seldom optimiim for any specific operating 
regime. Therefore, optimum performance does not reflect the 
fidelity of a motion system. Approximate duplication of 
such fundamental maneuvers as coordinated turn entries has 
been accomplished by matching aircraft and simulator accel- 
erations about all six degrees of freedom. Since a precise 
duplication is unattainable due to the simulator's physical 
constraints, only the human operator can properly evaluate 
the differences. 

Optimization of the control system for a new simu- 
lator is typically accomplished through numerous flights by 
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an experienced test pilot while engineers vary stability 
parameters. Many unpublished fidelity programs involve 
hundreds of manhours and expensive computer time in search 
of optimal control parameters. Often pilots need refresher 
flights in the aircraft to regain familiarity with its 
stability characteristics. One such study [Ref. 22] esti- 
mated that "eight hours of simulator time is sufficient to 
destroy a pilot's familiarity with the airplane to the point 
that he can no longer make valid qualitative judgements on 
flying qualities during closed loop tests." 

C. PRESENT PURPOSE 

Present D.O.D. emphasis is on new, higher fidelity 
flight simulators. A second avenue of approach is the 
improvement of present training devices having less than 
optimal flying qualities. Lack of fidelity has reduced the 
purpose of many devices to mere procedural trainers, while 
the motion and visual systems go unused. With new technology 
visual display systems, the shortcomings of the handling 
qualities and performance characteristics are magnified. 

As simulators assume a greater role in flight training and 
proficiency programs teaching advanced flying skills, the 
need for realistic duplication of the aircraft is apparent. 

The Naval Air Test Center (NATO) has been conducting a 
fidelity improvement program for the Navy's inventory of 
Operational Flight Trainers and Weapon System Trainers 
[Ref. 22] . A team of fleet experienced pilots, flight test 
engineers, simulator and computer experts have been 
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successful in identifying performance, stability, and con- 
trol parameters contributing to poor fidelity in all the 
subsystems. Although some hardware changes have been made, 
the majority of fidelity improvements have been accomplished 
by software reprogramming. 

"Aircraft manufacturers generally rely on wind 
tunnel data as a data base. Flight test data are not 
available because the first training device is usually 
delivered to the Fleet at about the same time that the 
first airplane comes off the production line. Since 
wind tunnel data represent at best only an estimate of 
airplane flight characteristics, the use of these data 
as a data base for a simulator results in poor flying 
qualities simulation. In addition, there has been no 
organized effort within the Navy or by any manufacturer 
to reprogram a simulator once flight test data are 
available . " 

"Recent advances in the practical applications of 
control and estimation theory have made it possible to 
extract stability derivatives and other parameters 
from flight test data to a degree and with accuracy 
not possible in the past." 

The degree of realism required in equipment, visual dis- 
plays, and motion systems is a function of the level of 
training being performed and the amount of learning transfer 
desired. For a specific training device, whatever cost- 
effective measures (up to the point of diminishing returns) 
increase its fidelity would be welcomed, thus improving the 
transfer of learning. 

The goal of this thesis is the development of a technique 
which would allow the simulator programmer (such as those of 
NATO's fidelity improvement program) to optimize the archi- 
tecture and choice of parameters for simulator software. 

Work on an ongoing research program by Dr. L. R. Young at 
the NASA Langley Research Center and the MIT Man-Vehicle 
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Laboratory on the integration of visual and motion cues 
provided many original ideas for this work. An experiment 
was performed at the NASA Ames Research Center in the Flight 
Simulation Laboratory on a six-degree-of-freedom research 
simulator. The experiment was designed to test the hypothe- 
sis that there was no difference between simulator motion 
cues presented to the pilots. Five pilots performed a total 
of 161 simulator "runs", providing subjective ratings of 
different degrees of fidelity. 

Using Ormsby's model of human dynamic orientation, a 
computer program was written to simulate the pilot's subjec- 
tive evaluation of fidelity in objective terms. Validation 
of the simulation was performed using the experimental re- 
sults as a data base. Another aim of this report is to 
confirm Dr. Ormsby's model to whatever extent possible. 
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II . METHOD 



A. EXPERIMENT 
1 . Design 

The pilot-subjects were asked to provide responses 
to different levels of motion fidelity. Fifteen sets of 
"wash-out" filter parameters were selected in an effort to 
produce discernably different motion cues during a flight 
task. For each "flight", subjects were asked to rate the 
rotational and linear acceleration cues relative to a 
rating scale and to provide a combined fidelity score. 

Thus, the dependent variables for the experiment were the 
three subjective ratings of the motion cues. The indepen- 
dent variables in the analysis were the motion cues and 
parameter sets. 

The experimental results were analyzed according to 
a three-way factorial analysis of variance. Two fixed fac- 
tors were motion (three levels) and parameter sets (15 
levels) . The third was subjects, a random factor (five 
levels) . A conceptual model of the experimental design is 
shown in Figure 2. 
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SUBJECTS S 3 
S. 




MOTION CUES 



Figure 2: 



Model of Experimental Design 



2 . Sub j ects 

Since the simulator was programmed as an F-15 (a 
combat fighter aircraft) , subjects for the study were chosen 
for their experience as pilots in high performance aircraft. 
Although the subjects obtained for the experiment were ex- 
perienced in military combat aircraft, only one was relative- 
ly current. The other four subjects were military trained 
airline pilots or flight engineers. Their pilot time exper- 
ience is shown in Table I. Each was a volunteer and was 
compensated for his participation. 
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SUBJECT HELP SINGLE -ENGINE MULTI-ENGINE HIGH PERFORMANCE JET 



1 


550 


2400 


3050 


2 


180 


860 


130 


3 600 


160 


350 


250 


4 


30 


- 


800 


5 


1600 


7500 


2150 


Table I; 


Subject 


Experience (pilot hours) 





3 . Stimuli and Apparatus 
a. Simulator 

The experiment was performed at the NASA Ames 
Research Center, in the Flight Simulation Laboratory. The 
simulator used in the experiment was the S.Ol All-Axis 
Motion Generator, manufactured by the Northrop Corporation. 
The six-degree-of-f reedom research simulator has independent 
drive mechanisms with a single-seat cockpit in the center of 
its operating envelope. The machine's linear position limits 
are ±9 feet; the rotational limits are ±45 degrees from its 
neutral position. The simulator was controlled by an EAI 
8400 computer. The computer programming included "soft" 
position, velocity, and acceleration limiting networks which 
prevent the simulator from reaching the "hard" mechanical 
limits. A more in-depth study of the machine's operating 
characteristics and mathematical development of its "drive 
logic" is contained in Reference 15. 

As with most simulators, the cockpit was en- 
closed preventing visual reference to outside of the cab. 
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The cockpit configuration was that of a typical single-seat 
aircraft, controlled by a stick, rudder pedals, and power 
lever . Since visual cues from external displays and arti- 
ficial attitude instrxomentation tend to interact with 
vestibular cues in a way not fully understood, use of these 
displays was minimized. (The simulator was equipped with a 
terrain modelboard visual system, but the projection equip- 
ment was off for the duration of the experiment.) The main 
instrument panel was organized in a conventional "T” arrange- 
ment as shown in Figure 3 . Although the basic instrument 
array included an artificial horizon (attitude gyro) , it 
was felt that its absence would impose an unusually high 
workload level on the subject and tend to lengthen the learn- 
ing process inherent in the task. Also, the more familiar 
the cockpit environment was to the subject, the quicker 
would be his adaptation to the task, allowing maximum 
attention to the motion cues. 

An intercom system, with a "hot" microphone for 
the pilot, allowed two-way communication between the subject, 
system operators, and experimenter. Verbal responses by the 
pilots were recorded on a cassette tape recorder. The ex- 
perimenter monitored the progress of the task on a strip 
chart recorder, displaying altitude, airspeed, vertical 
speed, and heading. 
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1 - Clock with Sweephand 

2 - Airspeed Indicator 

3 - Attitude Gyro 

4 - Altimeter 

5 - RPM Gauge 



6 - Angle of Attack Indicator 

7 - Horizontal Situation 

Indicator 

8 - Vertical Speed Indicator 

9 - Turn Needle and Ball 



Figure 3: Main Instrument Panel 



The "wash-out" filter, which interfaces the 
mathematical aircraft model with the "drive" program, 
is composed of six separate digital filters, one for each 
degree of freedom or "channel". Each second order, high 
pass filter is described by a linear gain, a natural 
frequency, and a dam.ping ratio with a Laplace transform 
equation of the form: 
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Ks- 



+ w2 



where: K = linear gain 

s = Laplace complex variable 
^ = damping ratio 
Wj^ = natural frequency 
Through consultation with NASA personnel, 15 
sets of parameters (18 param.eters each) were chosen to 
produce noticeably different levels of realism. Values of 
the filter parameters are presented in Appendix B. The 
fidelity associated with each set could not be predeter- 
mined, except that the range over which each parameter 
varied was thought to be sufficient to produce noticeable 
differences in motion cues, 
b. Task 

Initially, the experiment was designed to ana- 
lyze two different tasks in order to exam.ine any significant 
difference in subject responses. Due to time constraints 
and equipment malfunctions, the second task, an ILS (instru- 
ment landing system) approach to a waveoff, was eliminated 
from the study. 

Design of the task reflected many considerations. 
It was recognized that the task should be: 

■ similar to those used in RDT&E and training programs 
so that the maximum, extrapolation of results could be 
made 

■ relatively familiar to the subjects and well within 
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their capabilities to lessen the time spent mastering 
it 

■ well within the capability of the simulator to mini- 
mize machine error and breakdowns 

■ long enough to expose the subject to a full variety 

of motion cues and allow him ample time to subjectively 
evaluate them. 

A three and one-half minute "basic airwork" task 
was chosen to simulate the low frequency flight maneuvers 
common to the terminal phase of flight operations (climb-out 
after takeoff, descent prior to approach and landing, etc.). 
Subjects were provided enough information to rapidly learn 
the task and to execute it consistently from trial to trial. 

As shown in Figure 4, two versions of the flight pattern 
were devised. The direction of the initial turn was random- 
ized. During the second replication of each parameter set, 
the turn was made in the opposite direction, 
c. Rating Scales 

Subjects were asked to subjectively judge the 
fidelity of the motion cues encountered during the task us- 
ing two rating scales. The pilots scored separately the 
linear (or translational) and the rotational motion cues on 
a graphic eight-point scale. The cues were also rated in 
comparison to a "standard" motion (presented using the 
established set of wash-out parameters employed on the simu- 
lator) . Thus, the subjects rated the motion fidelity on 
both an absolute and comparative basis. 
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Figure 4: Flight Task Pattern 
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Figure 5; Component Rating Scale 
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Figure 6: Comparative Rating Scale 



4. 



Procediire 



After the subjects read the instructions (see Appen- 
dix A) , further questions were answered. The intent and 
goals of the experiment were discussed freely with the sub- 
jects. They were allowed to become familiar with the simu- 
lator cab during one practice run before the trial runs 
began. Each was given a pilot kneeboard containing a dia- 
gram of the task, the two rating scales, and a list of topics 
from which relevant comments were volunteered following each 
run. 

The sequence of wash-out parameter sets presented to 
the subjects was randomized; a different sequence was pre- 
sented to each subject. The "standard” motion was presented 
once at the beginning of each session (twice at the start of 
a subject's first session) and every fourth run thereafter. 

A typical session of 12 simulator runs is shown in Figure 7 . 

^0 ^2 ^8 ^15 ^0 ^11 ^14 ^6 ^0 ^1 ^5 ^13 

where: Pq = "standard" parameter set 

Figure 7: Typical Sequence of Parameter Sets 

No reference was made to the parameter set number or their 

values dxiring the trial runs. 

At the end of each task, subjects rated the motion 
cues and provided subjective comments relevent to the spe- 
cific run via the intercom. The parameter set was changed 
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and the simulator was reinitialized for the next run. After 
each session, subjects were allowed to review their responses 
and discuss the session (comments were not recorded) . 

Factors which were uncontrollable were the number of 
sessions and runs per subject. Due to equipment malfunctions 
and scheduling constraints the sessions and runs varied 
widely as sho^m in Table II. 



SUBJECT 



NO. SESSIONS 



1 

2 

3 

4 

5 



3 

3 

4 
2 
2 



* Includes "Standard Motion" 



Table II: Sessions and Runs by 



NO. RUNS* 
24 
37 
42 
26 
31 

Runs 

Subjects 



B. SIMULATION 

1 . Model Assumptions 

The subjective process by which a pilot compares the 
motion of a simulator with that of an aircraft is, obviously, 
a complex psychological one. If a model of this human com- 
parative method were constructed, it could predict a pilot's 
response to the fidelity of a simulator's motion character- 
istics. The structure of the computer simulation described 
in this section is based on the assumption (as proposed by 
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Young, Ref. 13) that the pilot's evaluation is a function of 
the difference between: 

■ the simulator motion cues as perceived by the pilot 
and 

■ the perceived motion of the aircraft as experienced 
by the pilot. 

Assuming that Ormsby's vestibular model of human dy- 
namic orientation provides a reliable estimate of motion 
perception, it was then a matter of deriving motion histories 
of the aircraft and the simulator performing "identical" 
maneuvers . 

2. Design 

During the experiment described in Section A, the 
output of the computerized aircraft model (see Fig. 1) was 
recorded on magnetic tape during each run. The data con- 
sisted of accelerations for each degree of freedom recorded 
during each cycle of the discrete programming. The cycle 
time used in the experiment was 50 milliseconds, chosen so as 
to be compatible with the input sampling rate of the dynamic 
orientation program. (Sufficient data was also recorded to 
document the specific pilot and run number.) 

The flight acceleration data was used as an input to 
the S.Ol simulator's drive program, a copy of which is used 
in this simulation. The simulation is a deterministic model 
with no random variables. Figure 8 is a block diagram of the 
simulation. Since the dynamic response of the motion platform 
contributes to the motion cues presented to the pilot, but is 
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Figure 8: Simulation Block Diagram 




not a factor in this simulation, the lead compensation net- 
work in the program is bypassed. The result is the closest 
approximation available to the actual motion of the cab. 

The outputs from this program are linear acceleration com- 
mands and rotational position commands of the platform's 
drive mechanism. A transformation then converts these sig- 
nals back to vectors describing the cockpit motion. The 
rotational position commands are then differentiated to 
produce rotational velocity vectors. 

Inputs to the human dynamic orientation program are 
discrete time histories of specific force vectors (in units 
of "G's") and angular velocity vectors (radians/sec.) in the 
coordinate system of the pilot's head. It is sampled at dis- 
crete intervals of 100 milliseconds; thus the simulator and 
aircraft motion vectors are sampled at every other time 
frame. Although this represents a loss of input information, 
considerable recomputation of many model parameters would 
have been necessary to alter the input or output rates. The 
output of Ormsby's model consists of vectors in the coordi- 
nate system of the subject describing his perception of 
rotational velocity, linear acceleration, and attitude (per- 
ceived vertical or "down" direction). The model's output 
rate is one sample per second. 

Motion vectors of the aircraft are computed from the 
flight data after the rotational velocity is formed by inte- 
grating the rotational acceleration. These vectors are then 
used as inputs to the dynamic orientation program to produce 
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perceived aircraft motion. Notation used in the block 
diagram is defined as follows: 



= aircraft model linear and rotational acceleration 
vectors, respectively 



SF^ , SFg = specific force vectors for the aircraft and 
simulator, respectively 



~ rotational velocity vectors for the aircraft 
and simulator, respectively 

/S /N 

^ = perceived linear acceleration vectors for the 
aircraft and simulator, respectively 







= perceived rotational velocity vectors for the 
aircraft and simulator, respectively 

DOVTN^ = perceived attitude vectors for the air- 
craft and simulator, respectively. 
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III. RESULTS 



A. EXPERIMENTAL RESULTS 

1 . Subjective Responses 

Verbal responses recorded at the end of each task 
allowed the experimenter to understand such considerations 
as the learning curve, workload, etc ., associated with the 
task. The act of commenting on each run gave subjects the 
opportunity to critically review the flight and to adjust 
their scores as opinions were formulated. This section is 
a brief summary of pilot comments pertaining to the 
experimental flights. 

Differentiation between the rotational and transla- 
tional motion cues was discussed with the subjects as the 
experiment progressed in hopes that more reliable responses 
would result. Following his fifth run, one si±>ject attempted 
to further define his task: 

"I'm to separate pitch and roll from sideways, hori- 
zontal, or longitudinal axis motion as linear, whether 
it be throttle acceleration or just sideways gusts." 

Several pilots commented that they were more confi- 
dent of their responses during second and subsequent ses- 
sions as they familiarized themselves with the task and 
learned more about the characteristics of the motion to be 
rated. One pilot mentioned that it was difficult to dis- 
cern whether motion cues were more realistic or his ability 
to control the simulator had improved. During a subject's 
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second session he commented; 



....the oscillations I was in the other day; at 
first I thought I was experiencing pitch oscillations. 

What they are is the linear, up and down, that I was 
feeling before; so I'm getting used to the motion.... 
at least understanding it better.” 

All the pilots noted that the workload level varied 
significantly from run to run as the motion cues changed. 
Abnormal cues often caused subjects to react improperly; a 
response which presumably would not have occurred in a fixed- 
base simulator. One subject noted the difference between 
the motion and cockpit instrument indications; 

"I was getting a definite rotational pitch sensa- 
tion and frequently corrected with [forward] stick.... 
which I didn ' t want to do . ” 

All subjects noted vibrations induced by the motion 
base mechanism. Normally turbulence masks these anomalies to 
some extent. Since it was not used in the experiment, the 
mechanical limitations tended to put a ceiling on the subjec- 
tive motion ratings. One subject was allowed a practice run 
with a minimal level of turbulence at the beginning of his 
second session. He commented; 

"... turbulence ... simply detracts from the ability 
to differentiate. . .how the controls are functioning. 

I think I'm less able to determine how much of it is 
...turbulence or control response with the turbulence." 

2 . Correlation of Ratings 

Scores from the three rating scales represent discrete 
data possessing ordinal properties. It can be argued that 
each rating scale attempts to define levels separated by equal 
intervals; it is not immediately evident that each subject 
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perceived the scales in this manner. Thus, it cannot be 
assumed that the data comes from an interval scale, an 
assiamption for use of parametric statistics. As a measure 
of the reliability of the data, nonparametric correlations 
were calculated. Scores from the standard motion were ex- 
cluded from this and subsequent analysis since subjects were 
told which runs were the standard and instructed to use it 
as a comparative reference for the other 15 parameter sets. 
An average score was formed for each subject/motion/param- 
eter set combination analogous to each "cell" of the experi- 
mental design (see Fig. 2) . Spearman rank correlation coef- 
ficients between subjects were computed for the linear, 
rotational, and combined ratings separately. The correla- 
tions with respective significance levels are shown in 
Table III. 

SUBJECTS 







1 


2 


3 


4 


5 


s 


1 


1.0 


5240* 


.0731 


. 6810*** 


.1264 


u 














B 


2 




1.0 


.1435 


.5668** 


.1722 


J 














E 


3 






1.0 


.2263 


.5666** 


C 














T 


4 








1.0 


.0250 


S 
















5 










1.0 






Significance 


Levels ; 


*p<.05. 


**p<.01, *** 


p<.005 



Table Ill(a) ; Correlation Between Subjects for 

Linear Motion Ratings 
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SUBJECTS 
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2 


3 


4 


5 


s 

u 


1 1.0 


.4857* 


1244 


.6081** 


.4367 


B 

J 


2 


1.0 


2823 


.5957** 


.4157 


E 

C 


3 




1.0 


-.0703 


. 3963 


T 

S 


4 






1.0 


.2485 




5 








1.0 




Table Ill(b); 


Correlation Between Subjects 


for 






Rotational 


Motion 


Ratings 
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SUBJECTS 

3 


4 


5 














S 


1 1.0 


. 5449 * 


. 5065 * 


. 7907 *** 


.2456 


U 

B 


2 


1.0 


. 5321 * 


. 6147 ** 


.4700 


J 

E 


3 




1.0 


. 4481 * 


.4968 


C 

T 

S 


5 

Significance 


Levels : 


* P <. 05 , 


1.0 

** p <. 01 , *** 


.2564 

1.0 

p <. 005 



Table III(c): Correlation Between Subjects for 

Combined Motion Ratings 



Overall agreement between the five subjects was com- 
puted using the Kendall coefficient of concordance as shown 
in Table IV. 



47 



RATINGS 



Linear Rotational Combined 

.4483* .4835* .5914* 

Significance Level: *p<.005 

Table IV: Coefficient of Concordance Among Subjects 

for Motion Ratings 

The ability of subjects to repeat ratings given the 
same stimuli was investigated. Due to time constraints only 
41 subject/parameter set combinations were replicated. Cor- 
relations between the first and second observations for these 
combinations is shown in Table V. 

RATINGS 

Linear Rotational Combined 

.4882* .4234* .5559** 

Significance Levels: *p<.005, **p<. 0005 

Table V: Spearman Rank Correlation Between Replications 

3 . Analysis of Variance 

Analysis of Variance techniques were used to deter- 
mine if different motion cues significantly affected the 
ratings and if interactions existed which would help explain 
the weak between-subject correlations obtained. Although 
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the discrete rating scale data violates the assumptions of 
interval data and normality of the error distribution, an 
approximate ANOVA was computed. Since the rating scale con- 
sists of essentially 15 levels separated by intervals de- 
signed to be equal, the assumption of interval data is not 
entirely untenable. Madill [Ref. 17] reports that the dis- 
tribution of the error term in the regression model for 
discrete rating scale data is often binomial in nature. 

Thus, the variance is not independent of the mean as re- 
quired for ANOVA analysis. The transformation which attempts 
to equalize the error variance for the binomial distribution 
is the arcsinV^ transformation. A transformation of the 
form: 




where: R = raw rating scale score 

R' = transformed score 

was used. The effect is minimal in the range 3.0 to 6.0 and 
tends to expand the ends of the scale. The ANOVA table is 
shown in Table VI. Due to the large percentage of missing 
data, the three-way interaction sum of squares is included 
in the error sum of squares with its degrees of freedom. 
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SOURCE 


d.f . 


3.S. 


M.S. 


F RATIO 


Subjects (S) 


4 


83.91 


20.98 


44.24* 


Motion Cues (M) 


2 


1.08 


.54 


.33 


Parameter Sets (P) 


14 


119.99 


8.57 


4.86* 


S X M 


8 


13.13 


1.64 


3.46* 


S X P 


56 


98.71 


1.76 


3.72* 


M X P 


28 


13.25 


.47 


.99 


Error 


235 


111.91 


.48 




Total 


347 


441.98 







significance level; *p<. 001 

Table VI : ANOVA Table for the Experiment 

A nonparametric analysis of variance was performed 
to support the results of the approximate ANOVA above. The 
subject/parameter set averages were tested using the Friedman 
two-way analysis of variance by ranks [Ref. 18]. The three 
motion cues were investigated separately. Results showed 
that both the linear and rotational motion ratings were sig- 
nificantly different at the .01 level. The combined ratings 
were different at the .001 significance level. 

A Duncan multiple range test with unequal replications 
[Refs. 19 and 20] was performed on the means of the parameter 
sets. Results using the .05 level of significance are shown 
in Figure 9. Any two means underscored by the same line are 
not significantly different. 
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X-H- 

6.0 



Figure 9: Parameter Set Averages 

B. SIMULATION RESULTS 
1 . Fidelity Models 

The available signals at the simulation output con- 
sist of vectors estimating the pilot's perceived motion and 
attitude. Each human dynamic orientation program (see Fig. 8) 
computes perceived linear acceleration, angular velocity, and 
attitude vectors. Algorithms for comparing these time vary- 
ing signals were constructed based on the assumption that the 
greater the difference between corresponding vectors for the 
aircraft and simulator, the less the fidelity. The following 
vector variables are available as inputs to the fidelity 
function. (Notation is consistent with that used by Young in 
Ref. 12.) 

^(t) = perceived aircraft linear acceleration 

/N 

^(t) = perceived simulator linear acceleration 
^(t) = perceived aircraft angular velocity 
^(t) = perceived simulator angular velocity 
DOW^(t) = perceived aircraft attitude 
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DOWN^ (t) = perceived simulator attitude 



where : 



a = aircraft, s = simulator, indicates a perceived 

quantity 

(t) indicates a time 
varying function 



Three different fidelity functions were constructed 
to be compared with the experimental data. The first is an 
average percent error model as suggested by Young: 

"When the [motion] sensations are clearly supra- 
threshold, the most likely candidate is just percent 
error, the ratio of perceptual error to the correct 
quantity. When the model indications for 'correct' 
perceptions are subthreshold, it seems more logical 
to assess a large penalty for errors that are large 
compared to the threshold value . " 

For each sample output (one per second) the following "cost 

to realism" indices were computed. 



lA A(t)| 

TT^ i ^ 

max||^(t)!, 
lA $(t)| 

r,^ I /N S " ^w^^^ 

max[|w^(t)l, 

|A DO^(t)| = Cj^(t) 

where; A = the difference between aircraft and 
simulator vectors 

|x| = amplitude of vector x 

thr = threshold value 

Threshold values used in the simulation were those 
computed in Ormsby's vestibular modelling thesis [Ref. 11]; 



(2a) 

(2b) 

(2c) 
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♦ • I 



= .005 "G's" 

w,, = .0146 radians/sec 



The vector DOV7N is a constant amplitude unit vector; thus, 
the cost index is simply the magnitude of the difference 
vector. Averages for each of the three cost indices were 
computed over each simulator run. Though the task was de- 
signed to last three and one-half minutes, subjects often 
used over four minutes to complete it. Aircraft acceler- 
ation data from each run, up to four minutes, was used as 
input to the simulation. The following average cost values 
were computed: 



C 



A 




(ti) 



(3a) 



^w 

n 

Cn = 



n i=i D i 



(3b) 

(3c) 



where : n = nxamber of output time frames 

Several subjects noted that during many runs there 
existed unexplained and unrealistic motion excursions about 
the rotational axis, but more pronounced in the translational 
cues. A measure using the same cost indices which is more 
sensitive to these rapid anomalies (if caused by the program- 
ming and not the hardware) is a root mean square (RMS) fig- 
ure. The second fidelity function computes RMS cost values 

as follows: 
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(4a) 




(4b) 

(4c) 



It is not unreasonable to assume that pilot responses 
may simply be proportional to the amplitude of the error vec- 
tor- The third fidelity function computes the average of 
each error vector; the nxxmerators of equations 3a, 3b, and 



3c . 



A = 



w = 




(5a) 



(5b) 



i=l 

A total of 133 simulator runs were recorded success- 
fully; 97 of which used wash-out parameter sets other than 
the standard set. Minimum, maximum, and average cost indices 
computed by the three fidelity functions are shown in Table 



DOWN (tj^)l 



(5c) 



VII. 
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MINIMUM 


AVERAGE 


MAXIMUM 


1 

lu 


.9959 


1.0788 


1.4830 


Cw 


. 5776 


.7603 


. 9457 


C 

D 


.0626 


. 0898 


.1243 


Ca 


.9983 


1.1428 


2.0106 


Cw 


. 6750 


.8309 


1.0671 




.0715 


.1045 


.1451 


A 


.0403 


. 0612 


.0911 


W 


.0108 


.0179 


.0282 


D 


. 0626 


. 0898 


.1243 



Table VII: Cost Index Values Computed by the Simulation 

2 . Experimental Optimization 

Response of the simulation to changes in the wash-out 
parameters was investigated to determine its adaptability to 
experimental optimization methods such as an optimal gradient 
search. If the simulation is to be an effective research 
tool for determining optimal control system variables, then 
its response characteristics must be of a tractable form. 

Time did not allow for a full exploration of its response 
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"surface" by an optimization routine, nor were manual methods 
employed to search for a minimum "cost" solution. 

Linear and rotational wash-out filter gains were 
varied over the range of the parameters used in the experi- 
ment. The response function was arbitrarily chosen to be: 









( 6 ) 



The linear gains, were fixed at different levels 

while the rotational gains, K^, were varied together from 
0.0 to 0.8 at intervals of .05. The response curves obtained 
are shown in Fig. 10. The basic parameter set used during 
these simulation runs was . 



C. SIMULATION VALIDATION 

Validation of the simulation consisted of comparing the 
subjective fidelity ratings with predictions from the simu- 
lation. It was then a question of which ratings and cost 
indices to analyze. Since the strongest correlations be- 
tween subjects were for ratings of the combined motion cues, 
the combined ratings were compared to the following unweight- 
ed linear combinations of the cost indices: 



^A 


+ 


Cw 


+ 




(7a) 


C, 


+ 


c 






(7b) 


A 




w 




D 


A 


-h 


W 


+ 


D 


(7c) 



Spearman rank correlation coefficients were computed for 
each subject and cost index combination. (Fidelity ratings 
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were subtracted from a constant to yield scores propor- 
tional to the cost to fidelity. ) The results are shown in 
Table VIII. 



_ _N0. 

+ ^w + + ^w + A + W + D RUNS 





^1 


.4147* 


.5138** 


.2673 


14 


s 

u 


^2 


. 2893 


. 3445* 


-.0403 


21 


B 












J 

E 


S 3 


. 3802** 


. 3986*** 


-.1489 


26 


C 












T 

S 


S 4 


. 0542 


.2965 


-.0599 


13 




S 5 


. 3217* 


.3227* 


. 3953** 


23 



Significance levels: *p<0.1, **p<.05, ***p<.025 



Table VIII: Spearman Rank Correlations Between Subjects 

and Simulation Cost Indices 



Weighted sums of the cost indices were computed by mul- 
tiple linear regression techniques in order to improve the 
correlations. Using the assumptions for parametric methods 
given for the analysis of variance above, a regression model 
of the form: 



^c “ ^l^w ^2^A ^4 

where: R^I, = 8 - R^, 

R^ = combined motion rating 
a = regression coefficients 



(S) 
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^liOOOOOB: 
EOOOOB5S 
Till a' 






was used to compute regression coefficients for each subject 



and for all subjects combined. Results are shown in Table IX. 



Regression Coefficients 
a. 



Sig. @ 



S 

U 

B 

J 

E 

C 

T 

S 



5 

All 
Subj ects 



-.530 


2.697* 


5.676 


-.874 


.5660 


.2571 


2.819 


1.845 


-17.628* 


. 547 


.4926 


.1826 


3.578 


2.030 


1.811 


-2.926 


.3671 


.3541 


-6.540 


.449 


-7.916 


9.312 


.2949 


.8346 


8.792* 


.780 


9.096 


-4.999 


.4589 


.2030 


-1.336 


1.513** 


-3.070 


2.653 


.2442 


.1245 



(§ Significance level associated with the ANOVA F ratio for the 
regression sum of squares 

Significance levels; *p<0.1, **p<0.05 



Table IX: Regression Coefficients for Combined Ratings 

versus RMS Cost Indices 
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IV. DISCUSSION 



Results of the experiment tend to reject the hypothesis 
that similar motion cues were presented to the subjects. As 
shown in Figure 9, a cluster of nine parameter sets exists 
toward the center of the range of means. Parameter sets one 
and three appear significantly better than the others, while 
sets 9, 11, 14, and 15 are significantly inferior to the 
others. The fact that the center group of parameter sets 
produced indistinguishable motion cues is a probable explana- 
tion for the weak correlations between subjects. Parametric 
Pearson product-moment correlations were computed yielding 
similar results at comparable significance levels. The 
highly significant coefficients of concordance also indicate 
that ordinal agreement between subjects was less than antic- 
ipated. The overall higher correlations for the combined 
motion ratings infer that subjects were capable of yielding 
more reliable scores after rating both linear and rotational 
motion cues on an absolute scale, followed by a pairwise 
comparison with the standard motion. The higher correlations 
could also be explained by the fact that several subjects 
said their responses resulted from an overall impression of 
the motion cues rather than a subjective average of the 
linear, rotational, and comparative scores. Two subjects 
rated the combined motion cues higher than either linear or 
rotational scores for several runs. 
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Comments by the pilots indicate that a large percentage 
of the variance of the pilot ratings is likely due to un- 
realistic dynamics of the motion drive mechanism. High 
frequency vibrations in the system were not as distracting 
and did not seem to lessen fidelity as much as the low fre- 
quency, larger amplitude anomalies experienced by the sub- 
jects. Whether these were caused by the software or 
hardware is not known. 

Cost indices computed by the simulation show large dif- 
ferences between outputs of the vestibular models. Errors 
ranged from 58 to 148 percent of the perceived aircraft 
motion vectors as shown in Table VII. Whether this typifies 
the capabilities of the S.Ol simulator or is largely a result 
of programming limitations and inaccuracies is not readily 
apparent. A significant amount of error is most likely due 
to roundoff, precision limits, and the discrete nature of the 
program. Motion histories of the aircraft and simulator have 
cycle times of .05 seconds, while the output rate of each 
vestibular model is one per second. This represents a con- 
siderable loss of information. If the output sampling rate 
equaled the input rate, it may have been more responsive to 
the rapid motion anomalies noted by the subjects. A quick 
excursion lasting less than a second may not have been re- 
flected in the output of the human dynamic orientation model. 

The aircraft model used during the experiment was that of 
the F-15 fighter. Ideally, subjects with recent F-15 experi- 
ence should have been chosen, but due to time and availability 
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a- 



fV 



I 



constraints selection of the pilots was less than optimal. 

It was felt that this factor contributed greatly to the lack 
of higher correlations of the ratings between subjects. 

This conclusion was supported by the weak correlations be- 
tween the first and second replications of the parameter 
sets. This also implies an inability of the subjects to 
repeat the task as well as inconsistent responses between 
subjects . 

Response of the simulation to changes in wash-out filter 
gains indicated a definite tendency for convexity in the 
neighborhood of minimum output. This can be explained in- 
tuitively. VThen the gains are relatively low, the lack of 
motion cues produces a large discrepancy between perceived 
aircraft and perceived simulator motions. When the gains 
are high, cues become overly rapid, exagerated, and fre- 
quently cause the simulator to reach its "soft" programming 
limits. When reached, the limits force the linear and ro- 
tational cues to become uncoordinated, reducing fidelity. 

The simulation output, defined by equation 6, varied over a 
small range as the rotational gains were varied as shown in 
Figure 10. For each level of the linear gains, the output 
varied only 2 to 4 percent. Parameter sets which produced 
significantly better motion cues, and P^, were comprised 
of filter gains set at 0.3. The tendency for the simulation 
output to minimize at approximately 0.3 indicates a capacity 
to predict subjective preferences for motion cues produced 
by these parameter sets. The simulation was also run with 
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the linear gain, , set at 0. Response was less than that 
for set to 0.1 throughout the range of K^. This implies 
that optimal motion cues include very little linear acceler- 
ation for the specific flight used as input to the simulation. 

Comparison between the simulation output and subject 
responses is summarized as follows; 

S Positive correlations were obtained for the first two 
fidelity models, equations 7a and 7b (see Table VIII). 

■ The third function, equation 7c, showed little correla- 
tion with the exception of subject five. 

■ The second fidelity function consistently produced 
higher and more significant correlations infer ing that 
the RMS cost indices are a better method for simulating 
subject responses. 

■ Linear regression using the RMS cost indices as inde- 
pendent variables and the combined motion cue ratings 
as dependent variables (equation 8) showed little con- 
sistency of regression coefficients among the five 
subjects (see Table IX) . 

■ Linear regression analysis yielded multiple correla- 
tions not appreciably better than the correlations 
for the unweighted cost index combination. 

■ Very little correlation was found between ratings of 
the linear and rotational motion cues and the corres- 
ponding cost indices, such as the rotational ratings 
and 
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Nine of the 15 wash-out filter parameter sets were 
grouped in the center of the range of mean responses as 
shown in Figure 9. Since these did not represent signifi- 
cantly different means, Spearman rank correlations were com- 
puted between the subjects' combined scores and the second 
fidelity function, equation 7b, for the remaining six param- 
eter sets. The results are shown in Table X. 



SUBJECTS 




Correlation .9710** .0748 .6545* .5000 .7660* 

No. of Runs 6 7 10 3 8 

Significance levels: *p<.05, **p<.01 

Table X: Spearman Rank Correlations for the Remaining 

Six Parameter Sets 

With the exception of subject two, the subjects showed a 
marked increase in correlation with the simulation results. 
Although the sample sizes for this post hoc analysis were 
small, three of the five correlations were significant. It 
was noted that the correlations compared closely, rankwise, 
to the flight experience each pilot had in high performance 
jet aircraft as indicated in Table I. 
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VI . RECOMMENDATIONS 



Recommendations for further research include suggestions 
for improvement of the experimental data base and the com- 
puter simulation. Better reliability and consistency in 
subject responses and a more accurate and responsive simu- 
lation would improve the prediction capability of the 
simulation. 

Further studies should employ subjects proficient in the 
simulated aircraft. Since the results appeared to be strong- 
ly affected by the subjects' ability to repeat the task, use 
of experienced test pilots familiar with flight dynamics 
would help create a more reliable data base. Future re- 
search should include preliminary studies to determine other 
control parameters which when varied would create motion 
cues more significantly different. Discrete rating scale 
data is difficult to analyze with classic analysis of vari- 
ance techniques, therefore, methods for obtaining pilot re- 
sponses from a continuous, interval scale should be pursued. 

Time did not allow for elimination of the inefficiencies 
of the simulation. Recomputation of the human dynamic orien- 
tation program's input and output model parameters would most 
likely improve the responsiveness of the model. Converting 
the program to double precision would reduce the inaccuracies, 
but double the computer memory requirements. The core re- 
quired as listed in Appendix C and as executed on the IBM 
360/67 Operating System was 458K. The entire flight 
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acceleration data set was read in at once. Other time his- 
tories, such as the simulation motion vectors, were computed 
creating large core requirements. If each time frame were 
processed through the simulation before the next was read in, 
the core requirements would be reduced by 75 percent. 

It was assumed that an accurate facsimile of the simu- 
lator motion was generated by the drive program. Since the 
lead compensation network was bypassed in the simulation, 
follow-on studies should compare the output of the drive 
program with actual simulator motion histories. This is a 
possible source of considerable error. 

The motion vectors for the simulator and aircraft could 
have been compared directly without use of the human dynamic 
orientation program. A similar analysis of the difference 
between these vectors would provide a way of validating 
Ormby ' s model. If the method incorporating the model cor- 
related appreciably higher with subject ratings, a reasonable 
validation of the human dynamic orientation would be achieved. 
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APPENDIX A 



INSTRUCTIONS TO SUBJECTS 

A. PURPOSE OF THE EXPERIMENTS 

The intent of this study is to analyze pilot response to 
varying degrees of simulator realism or fidelity. Subjective 
ratings by commercial pilots like yourself will help con- 
struct a mathematical model of the human process of compar- 
ing simulator motion characteristics with those of the 
simulated aircraft. Ultimately, our goals are to validate a 
model of human dynamic orientation and to optimize the drive 
programs of aircraft simulators. 

B. TASKS 

You will be asked to repeat two familiar tasks, each 
lasting from three to three and one half minutes. The first 
will be a basic airwork pattern, designed to employ the 
maneuvers common to operations in the terminal phase of 
flight. The second task is a full ILS approach to a wave- 
off. The basic airwork task will be repeated approximately 
fifty times, while the drive program parameters are varied 
to produce different degrees of fidelity. At the end of 
each task you will be asked to rate the characteristics of 
the motion using three rating scales discussed below, as 
well as provide your subjective opinion of the motion, task, 
etc. During subsequent sessions, you will rate the motion 
following the ILS approach task, repeated roughly the same 
number of times. 
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Each session will consist of approximately 12-15 flights 
and last no longer than 80 minutes. Detailed diagrams of 
each task will be available during the runs. Additional 
details will be provided before the sessions begin. 

The simulated aircraft in this study is the F-15, a high 
performance, fighter aircraft. Although you probably have 
had no experience in the F-15 itself, your background in 
high performance type aircraft was a major consideration in 
your being selected as a pilot for this study. 



C. USE OF THE RATING SCALES - OBJECTIVE RESPONSES 

Each trial run will employ a different set of drive pro- 
gram parameters. A standard set of parameters will produce 
a "standard motion" which will be presented during the first 
run and every fourth run thereafter. Repetition of the 
standard motion will reinforce your perception of it, allow- 
ing you to make comparative evaluations of the trial runs. 

A typical sequence of 12 runs per session will be presented 
as follows: 



Trial No. 



0 



□ 



10 



11 12 



where 



5 indicates that the fifth run of the session will 



be the "standard motion." 

1 . Component Rating Scale 

Evaluate the linear (translational) and rotational 
components of the motion following each simulator run using 
the component rating scale. Rate the motion characteristics 
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based on the descriptive phrases provided at each level of 
the scale. These subjective descriptions should be used 
merely as an aid in your evaluations. Responses should 
reflect consideration of: 

a. Response to Control Inputs - simulator 
reaction to movement or pressure on flight controls. 

b. Motion Deficiencies - motion deviations from 
experienced or anticipated aircraft motion, not a result of 
control inputs (e.g. response to turbulence, oscillations, 
etc . ) . 

c. Other - any other discrepancies in the motion 
dynamics which you detect and consider significant. 

Your responses should be formulated in reference to 
the absolute scales as well as in comparison to your rating 
of the standard motion. 

2 . Comparative Rating Scale 

Rate the trial runs in comparison the the standard 
motion using the comparative rating scale. Brief descrip- 
tions of each level of the scale should aid your evaluation. 

3 . Combined Rating Scale 

Using the descriptive phrases on the component rating 
scale, rate the trial run on an overall or combined basis. 
Reference should be made to your scores given on both the 
component and comparative rating scales. 

4 . General 

a. If you cannot choose between two adjacent levels 
on either of the scales, a "half" level score is allowed. 
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provided careful consideration is given to both adjacent 
levels. 

b. Ensure that your rating of the "standard motion" 
is recorded on the scale and referred to while rating trial 
runs . 

c. Basic performance measures will be recorded 
only during the ILS approaches. Although your performance 
will be compared to that of other pilots, our interest lies 
in how your performance varies with changes in fidelity. 

d. Try to ignore changes in your performance on 
the task as each session progresses. 

e. Your evaluations should reflect your assessment 
of all the motion cues detected. 

D. SUBJECTIVE COMMENTS 

After rating each trial run, brief comments on the items 
listed below are encouraged. It is not necessary to remark 
on each issue, but only on those you feel significant or 
particularly relevant to the trial run. 

1. Your ability to rate, separately, the linear and 
rotational components of the motion. 

2. Subjective description of the motion characteristics. 

3. Difficulties encountered with the task as it is 
affected by the particular motion characteristics. 

4. Impressions of the evaluation task. 

5. Duration of the task versus your ability to evaluate 
the motion. 
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6. Physical discomfort (fatigue, nausea, etc.). 

7. Workload level in terms of the attention level 
required to perform the task. 
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APPENDIX B 



WASH-OUT FILTER PARAMETER SETS 

Each wash-out filter (one for each degree of freedom) 
is described by three parameters. The values of the param- 
eters used in the experiment are shown in Table Bl. Where 
only one figure is shown, that value was used for all three 
channels. Natural frequencies for the translational chan- 
nels are in the order of roll, pitch, and yaw. The 
parameters are defined as follows: 

- Linear gains for the translational channels 

- Linear gains for the rotational channels 

w^ - Natural frequencies for the translational channels 

w - Natural frequencies for the rotational channels 
r 

- Damping factors for the translational channels 

- Damping factors for the rotational channels 

The "standard motion" parameter set is denoted by Pq. 
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Parameter 



Set No. 




K 

r 






"^r 






^0 


0.5 


0.5 


.4 


.4 1.0 


1.2 


. 707 


. 707 


^1 


0.3 


0. 3 


.2 


.2 .2 


0.7 


.707 


.707 


^2 


0.1 


0.1 


.4 


.4 1.0 


1.2 


.707 


.707 


^3 


0.3 


0.3 


.2 


.2 .7 


1.6 


. 707 


. 707 


^4 


1.0 


0.1 


. 7 


.7 1.5 


0.7 


.707 


. 707 


P5 


1.0 


0.1 


. 7 


.7 1.5 


1.2 


.707 


.707 


^6 


1.0 


0.1 


.7 


.7 1.5 


1.6 


.250 


.707 


"7 


0.7 


0.3 


. 7 


.7 1.5 


1.2 


.707 


. 707 


^8 


0.1 


1.0 


.7 


.7 1.5 


1.6 


.707 


.707 


^9 


0.3 


0.7 


.2 


.2 .7 


0.7 


.707 


.707 


0 

1 — 1 
0 ^ 


0.3 


0.7 


.2 


.2 .7 


1.2 


.707 


.707 


1 — 1 
1 — 1 
0 ^ 


0.1 


1.0 


. 7 


.7 1.5 


0.7 


.707 


.707 


^12 


0.1 


1.0 


.7 


.7 1.5 


1.2 


.707 


1.5 


"13 


0.3 


0.7 


.7 


.7 1.5 


1.6 


.707 


.707 


^14 


0.7 


0.7 


. 7 


.7 1.5 


0.7 


. 707 


.707 


Pl 5 


1.0 


1.0 


.4 


.4 1.0 


1.2 


.707 


.707 


Table 


Bl: Wash 


-out 


Parameter Sets 


Used in 


the Experiment 
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APPENDIX C 



SIMULATION SOURCE PROGRAM LISTING 

The simulation, written in FORTRAN IV H, consists of sub- 
routines which perform the functions shown in the block dia- 
gram in Figure 9. The program listing in this appendix is in 
the configuration used to perform the experimental optimiza- 
tion described in section III.B.2. The main subroutine 
(called SIM) controls the other subroutines and prints the 
wash-out parameter set and messages displaying the progress 
of the program execution. A main program (not listed) reads 
in the values for the wash-out parameters, varies them, and 
iteratively calls the main subroutine. The flight acceler- 
ation time histories are read in by subroutines READ. The 
model parameters used by the human dynamic orientation pro- 
gram (subroutine HDO) are those obtained with the program 
listing contained in Ref. 21. The simulation was executed 
on an IBM 360/67 Operating System using single precision, 
floating point variables. 

Subroutines MCPY and MINV are IBM system library sub- 
routines called from subroutine COTRN to perform matrix 
operations. MCPY copies one array variable into another. 

MINV inverts a square matrix. 
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:^s:^c;4s:^s:4c3{c3}e :^c 3ic:}c :4c jje :$c :^C3}c :Jca;c jje:}: a^csjc :*c24c3*c3jc :^c:4«:^c3;c5^:Jc :{c :4c3jc :4c ^sjcsjcjjc :{c :^c:4: :4c 3|c 



:4c :4c :Jc 3 jt: 4 «, 3 |c 34 c :41:4c :4c :4c 



MAIN SUBROUTINE 



* *nr^ ' 



c:4c:4c34s 3}ti 



C3{c:4c2is:<< 



:;c :4c :4:34 s :4c :4c ^!c :4c 3}c 



-'f* *1 



c 

c 

c 

c 

c 



c 



c 

c 

c 

c 

IOC 

c 



c 

c 

c 



c 

109 

c 

c 

c 

131 



c 

c 

c 



PROGRAM TO EVALUATE THE SOFTWARE PARAMETERS OF thF S.Ol 
“^CTICN SIMULATOR USING THE ORMSBY HUMAN DYNAMIC 
ORIENTATION MODEL 

SUBROUTINE S I M { C A Dt CA VG ,CRM S , K ) 

CGMMON/XFLOAT/TMT (3,3) , TMR ( 3 , 3 ) ,BTACC(3) ,BRACC{3) 

=«', TDDS( 3) ,RSFX( 3 ) 

CCMMON/ IFIXED/IMODE, lOUMFU, I FORD, ICHG 

COMMON/ WSHO UT / AKR ( 3) ,AK1( 3) ,OMGR( 3),0MGT (3 ) ,ZETR(3) 

*, ZETT(3) 

CCMMCN/FL IGHT/XA( 4800,3) , XR ( 48 00 , 3 ) , P , C , I P AR 

DIMENSION ACC (3 ) , ROT ( 3 ) , ANG ( 4800 , 3 ) 

DIMENSION SFS( 4800,3) ,WS(4800,3) , WA(4800 ,3 ) 

DATA DT/.05/,N/4800/ 

PRINT THE WASHOUT FILTER PARAMETERS 

WRITE(6,100) ^ ^ . 

FORMAT! IHO, 'WASHOUT PARAMETERS READ FROM DATA FILE:') 

WRITE( 6 , 120 ) ( AKR( I ) ,1 = 1,3) 

WRIT£(6,12l MAKl ( I), 1 = 1, 3) 

WRITE! 6, 122)!OMGR!I ) ,1 = 1,3) 

-^RITE!6,123)(GMGT! I ), 1 = 1,3) 

.nPITE!6,124)!ZETR!I) ,1=1,3) 
WRITE!6,125)(ZETT!I),I=1,3) 

READ THE FLIGHT DATA (AIRCRAFT MODEL OUTPUT) 

I F (K.GT .1 ) GO TO 131 
CALL READ(NSP) 

FCRMAT!* • ^°FLI3HT DATA READ FROM TAPE FILE') 

INITIALIZE THE INPUT ACCELERATIONS TO ZERO 

BTACC ! 1 )=0. 

ETACC! 2 )=0. 

ETACC! 3)=-32.2 
BRACC! 1 ) =0. 

ERACC! 2)=0. 

ERACC!3 )=0. 

SET THE "FLAGS" FOR SUBROUTINE SOIMOT 



I 0UMFU=0 

IFCRD=0 

ICFG=1 

C SET TEE MODE TO "INITIALIZE" AND CALL SUBROUTINE SOlMCT 
C 

I MCOE=- 1 

CALL SOIMOT(DT) 

C 

101 FCRMAt!' ^^SUBROUTINE SOIMOT INITIALIZED') 

C SET TEE MODE TO "RUN" 

^ I MODE = 1 

C SET THE INPUT ACCELERATION SAMPLE AND CALL SUBROUTINE 
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W 4<- -ii- -M- 



c 

c 



c 

c 

c 

c 



c 

ZQ 

C 

102 

C 

C 

C 

C 

C 



103 

C 

C 

C 

c 



104 

c 

c 

c 

c 



105 

C 

C 

c 

c 



107 
C 

C 

C 

132 

C 

108 
C 



SOIMCT TO CREATE THE DRIVE COMMANDS 

DC 20 1=1, NSP 

BTACC( n=XA(I, 1 ) 

BTACC(2 )=XA (1,2) 
eTACC(3)=XA(I,3) 

8RACC( 1 )=XR (1,1) 

ERACC(2 )=XR (1,2) 

8RACC(3)=XR(I ,3) 

CALL SOIMOT(DT) 

TRANSFORM THE DRIVE SIGNALS BACK INTO CAd (PILCT) COORDI- 
NATES 

CALL COTRN ( TDDS ,TMT, 1 , ACC ) 

CALL C0TRN( RSFX,TMR,-1 ,ROT> 

SFS(I,1 )=ACC(1) 

SFS( I ,2 )=ACC( 2) 

SFS( 1,3 )=ACC( 3)-32.2 
ANG( I ,1 )=ROT( 1 ) 

ANGU , 2 ) = ROT( 2) 

AN6( l,3)=RGT( 3) 

CCNT INUE 

'^RITE( 6,102 ) 

FCRMAT(' SIMULATOR DRIVE SIGNALS CREATED') 

DIFFERENTIATE THE ROTATIONAL POSITION COMMANDS TO OBTAIN 
ROTATIONAL VELOCITY FOR USE BY THE HUMAN DYNAMIC 
ORIENTATION (HOO) SUBROUTINE 

CALL DIFF(DT,ANG, WS,NSP,N) 

WRITE(6,103 ) 

FCRMAT( • ROTATIONAL DRIVE SIGNALS DIFFERENTIATED') 

CALL SUBROUTINE HOO TO COMPUTE THE PERCEIVED SIMULATOR 
MOTION VECTORS 

CALL HOC( OT ,NSP ,SFS,WS,9,N) 

WRITE (6, 104) 

FORMAT! • PERCEIVED SIMULATOR MOTION VECTORS CREATED') 

INTEGRATE THE ROTATIONAL FLIGHT ACCELERATION TO BE USED 
BY SUBROUTINE HOO 

IF(K.GT.l) GO TO 132 
CALL IMT(OT,XR,WA,NSP,N) 

WRITE(6,105) 

FORMAT!' ROTATIONAL FLIGHT ACCELERATION INTEGRATED’) 

CALL SUBROUTINE HOO TO COMPUTE THE PERCEIVED AIRCRAFT 
MOTION VECTORS 

CALL HOC(OT,NSP,XA,WA,S,N) 

WR ITE (6 , 107 ) 

FORMAT!' PERCEIVED FLIGHT MOTION VECTORS CREATED') 
OPERATE THE COST ANALYSIS ROUTINE 
CALL CCST(CAD,CAVG,CRMS ) 



WRITE(6,108) 

FCRMAT(1H0,’ COST ROUTINE COMPLETE') 
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C FCRMAT STATEMENTS 
C 

120 FORMAT! /5X, ‘AKR( 1 ,2,3) 

121 FGRMAT(5X,»AK1(1,2,3) 

122 FORMAT! 5X, ’OMGR( 1 ,2,3) 

123 FORMAT! 5X, ‘OMOT! 1 ,2,3) 

124 FORMAT! 5X, • ZETR! 1 ,2,3) 

125 FORMAT! 5X, ’ ZETT! I ,2,3) 
C 

RETURN 

END 



- I 



— , 



• ,3F6.3) 
, 3F6.3 ) 

, 3F6.3) 
,3F6.3) 

, 3F6.3 ) 

, 3F6. 3/ ) 



W ^ ^ ^ ^ % 

^ -nr -!• -T ^ 



SUBROUTINES 



: 3^ :4s 3jt : 



s^: 4s3j: 4 s:4::4s3*s:^:4s :4s3^3ic5js-3;s3jc:is^;s3»c55s:4s:4s:4: 



c 

c 

c 

c 

c 



c 



c 



SUBROUTINE TO INTEGRATE THE TIME SERIES CONSISTING OF THE 
ROTATIONAL FLIGHT ACCELERATION (AIRCRAFT MODEL OUTPUT) 



SUBROUTINE INT ( H , Y , Z ,NSP , N ) 
DIMENSION Y(N, 3) ,Z(N,3) 



HT=.3333333*H 
00 7 J= 1 , 3 
SUM1=Y( 2 ,J ) +Y( 2 , J ) 

SLM1=SUM1+SUMI 

SUMUHT’f^! Y! 1 , J ) +S UM l+Y ( 3 , J ) ) 

ALXl=Y( 4, J ) +Y ( 4, J ) 

ALX 1=AUX1+AUX1 

ALX1=SUM1+HT«!Y(3, J)+AUX1+Y( 5, J) ) 

AUX2=HT*! Y( l,J)+3 . 375* ! Y ( 2 , J ) +Y ! 5 , J ) ) +2 . 62 5 *! Y ( 3, J ) + 
*Y(4,J) ) +Y! 6,J) ) 

SUM2=Y( 5,J )+Y(5, J ) 

SUM2=SUM2+SUM2 

SLM2=AUX2-HT*(Y(4,J)+SUM2+Y (6, J) ) 

Z! 1 , J)=0. 

AUX=Y( 3 , J) +Y! 3, J ) 

ALX=AUX+AUX 

Z( 2, J) = SUM2-HT*(Y (2, J) +AUX+Y(4, J) ) 

Z! 3,J) = SUM1 
Z(4,J)=SUM2 



DO 4 I=7,NSP,2 

SLM1=AUX1 

SUM2=AUX2 

ALX1=Y! I-l, J)+Y( I-1,J) 

AUX1=AUX1+AUX1 , 

AL'X1=SUM1+HT«! Y! I -2 , J ) +AUX1 +Y ( I , J ) ) 

Z ( 1-2, J )=SUM1 
IF! I-NSP) 3 ,6,6 
AUX2=Y! I,J)+Y! I,J) 

AUX2=AUX2+AUX2 

ALX2=SUM2+HT*!Y( I - I , J ) + AUX 2 +Y ! I+l, J) ) 

Z! I-l, J )=SUM2 

CCNT INUE 

Z(NSP-1 , J)=AUX1 

Z(NSP, J)=AUX2 

GO TO 7 

Z!NSP-1 ,J) =SUM2 
Z (NSP, J )=AUX1 
CONTINUE 



RETURN 

END 
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OO O O O O o OOOOOOO 



C SLSRCUTINE TO PRODUCE THE SIMULATOR DRIVE SIGNALS FROM 
C TFE FLIGHT DATA OBTAINED FROM THE MAIN PROGRAM 
C 

SUBROUTINE SOIMOT(DT) 

MOTION DRIVE SYSTEM FOR S.Ol (SIX-OEGREE) 

PREPARED BY J.OOUVILLIER (NASA) 3 SEPT 1975 

PROGRAMMED BY O.ASTILL (CSC) 12 SEPT 1975 

MODIFIED BY DATA ANALYSIS GROUP (CSC) DEC 1976 

CCMMQN/XFLOAT/TMT ( 3, 3 ) ,TMR{ 3,3) ,BTACC(3) ,BRACC(3) 
T00S(3) ,RSFX(3 ) 

CCMMON/IFIXED/IMODE, IDUMFU, IFORD, ICHG 
CCMMQN/ WSHOUT/AKR(3 ) , AK 1 ( 3 ) , OMGR ( 3 ) ,OMGT (3 ) , ZETR( 3 ) 

*, ZETT(3) 

CCMMON/MOTION/QMGRL( 3) ,OMGRQ( 3) ,OMGRE( 3) ,OMGTE( 3) 

1, ZETAQ(3),ZETRL(3),ZDLF(3) , GNR S ( 3 ) , AKRL ( 3 ) ,SOMGR(3) 

2, Z0MGR(3),0MGQSQ(3),T0MGRE(3) ,OMGR£2(3) ,OMGT2(3) 

3, TOMZT(3) ,F0MGTE(3) ,SOMGT2(3) ,FOMGT3(3) ,0MGTE4(3) 

4, CM.GTQ(3) ,GLEAD(3 ) 

DIMENSION RPOSL (3 ),TPOSL(3) , TMS ( 3 ) , RAC ( 3 ) , T AC ( 3 ) 

1 , FNOT( 3 ) ,T0CS I( 3) ,TDCS 11(3), FNCHI ( 3 ) , FNC H 1 1 ( 3 ) 

2, RHFI (3),RHFII(3j ,RDHI ( 3 ) , RCH ( 3 ) , R OH ( 3 ) , FN A ( 3 ) 

3, RHFF3( 3) ,R0L I (3) ,RLF( 3) ,RL( 3 ) ,RSF(3 ) ,RDM( 3 ),RDCH( 3) 

4, FCF(3) ,RHF(3) ,RS (3) ,RIC(3) , T I C ( 3 ) , RSD ( 3 ) , FNCL ( 3 ) 

5, TDSC(3 ) , FCHK 3),FKN(3) ,FNAS(3) , FNCHI 3) ,FNChTP(3) 

6 , TDDH(3 ) ,TOOT (3 ) , ADDL( 3) , TDS ( 3 ) , TS IL ( 3 ) , TOCSP ( 3 ) 

7, TODA(3) ,TDDG(3) ,TSO(3) ,AFU(6) , ADSD ( 6 ) , K FL A G ( 3 ) 

8, FOLUPS(6) ,GNH(3),AIC(6) 

equivalence ( FOLUPS( I ) , AF'J( 1 ) ) 
equivalence ( ADS0( 1),TSD( 1) ) 

DATA G/32. 2/, RG/. 031056/ 

CATA TODG/0 ., 0. ,-32.2/ 

DATA FNOT/3’^1./ 

CATA TIC/3=i'0./,RIC/3=i=0./ 

DATA AODL/5.6,6 .8,5.5/ 

DATA GNH/3*1.0/ 

DATA AKXL/1./ 

DATA AKYL/1./ 

DATA RP0SL/3*.5326/ 

DATA TPOSL/3*8./ 

DATA NORES/0/ 

FAST SINE AND COSINE APPROX. 

XCQS(X) =1 .-(X*X)/2. 

XSIN( X) =X-( X**3) /6. 



IF ( IMODE) 10 ,60,100 



10 DC 20 IN=1,3 
C 

KFLAG ( I N )=0 
FOLUPSI IN)=0. 
FGLUPS ( IN + 3 ) = 0. 
ACSD( IN )=0. 
ACSD( IN + 3)=C. 
RLF( IN)=0. 

TSC( IN) =0. 

TSIL( IN)=0. 

RCH( IN)=0 . 

RL( IN)=0. 

RDH I( IN )=0. 
RFFFBI IN)=0. 
DTP=0. 

C 

CPHSF=1 . 

CTHSF=1 . 

CPSSF= I . 



CLEAR LIMIT FLAGS, INITIALIZE 
FILTERS, AND ZERO DRIVES 
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c 

c 

c 

c 



2C 



2 4 
25 



c 

100 



51 

53 



c 

c 

c 

2CC 



SPhSF=0. 

STHSF=0 . 

SPSSF=0. 

CPHLF=1 . 

CTHLF = 1 . 

CPSLF=1 . 

SPHLF=0 . 

STHLF = 0 . 

SPSLF=0. 

TMS( IN) =. 5/A0DL( IN) 
CCNTINUE 



IF( ICHG.EQ. 1 ) GO TO 25 
IF(DTP-DT) 24»100t24 
ICHG=1 
DTP=0T 

CALL SOIMIC(DTP) 

ICHG=0 

COMPUTE FOLLOW UP IF NO CA3 MOTION 

DC 53 J=l,6 

IF( IDUMFU.EQ.O) GO TO 51 
AFU(J)=ADSO(J) 

FCLUPS(J)=AOSO( J ) 

GO TO 53 

AFU( J)=FOLUPS( J ) 

CGNT INUE 

FNATS=0 . 

DC 200 IR=1,3 

COMPUTE LCAD FACTORS 

FNA( IR)= 3TACC( IR )*RG 
FNAS( IR )=FNA( IR)=:‘FNA( IR) 

FKN( IR)=AK1 (IR)*FNA( IR) 

FNATS=FNATS+FNAS ( IR) 

ROTATIONAL WASHOUT 

RDM( IR) =AKP ( IR) -BRACC ( IR) 

RDCH( IR )=RDM( IR)-ZOMGR( IR)=i'RCH( IR)-SOMGR( IR)*RHFFB( IR) 
RCF(IR)=RCH(IR)+RDCH( IR)*0T 



UPDATE FILTERS IF ICHG=1 
OR DT CHANGED 



CGNT INUE 



TMR( 1 , 1 )=l . 

TMR( 1 , 2 ) =-CPHSF*SPSSF/CPSSF 
TMR( 1,3 )=SPHSF^SPSSF/CPSSF 
TMR(2,1 )=0. 

TMR( 2, 2 ) =CPHSF/CPSSF 
TMR(2,3 )=-SPHSF/CPSSF 
TMR(3, 1 ) =0, 

TMR( 3, 2 )=SPHSF 
TMR( 3 , 3 )=CPHSF 



rotational (CAB TO GIM0AL) 
TRANSFORMATION 
MATRIX 



RCH(1) = RCH( 1 )+RCH( 2)*TMR( 1 , 2 ) +RC H ( 3 ) *TMR ( 1 ,3) 
RDH( 2) = RCH( 2)*TMR(2,2 ) +KCH ( 3 ) =?'T MR ( 2 , 3 ) 

ROH( 3) = RCH( 2)*TMR( 3,2)+RCH( 3 ) =!=TMR ( 3 , 3 ) 



FNAT=SQRT ( FNATS) 

I F( FNAT .EQ .0. ) FNAT=1.0 

FNMXL = FKN(1 )*AKXL/FMAT 
FNMYL=FKN(2 )*AKYL/FMAT 

FNMZL =-SQRT( I . - FNMXL^FNMXL-FNMYL^FNMYL ) 

LCW FREQ RESIDUAL TILT RATE 
RDLI(1)=(-FNMXL*( SP SL F«CT HL F ) / C P SLF - FNM YL *CPHLF*C THL F 
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c 

c 



c 



c 

c 



+FNMZL*SPHLF*CTHLF ) 

RDLI (2 ) = ( FNMXL*CTHLF/CPSLF +FNMY L*SPHLF ^S'STHLF 
+FNMZL^CPHLF*STHLF) 

ROLI (3 ) = ( - FNMXL*SPSLF^STHLF-FNMYL’i‘CPHLF-CPSLF*STHLF 
+FNMZL*SPHLF=*=CPSLF*STHLF) 

DC 250 IR=l,3 



RHFFB( IR)=RHFFB( I R ) +R0H ( I R ) =*=0T 
IF( IMOOE.LT.O) RDLI(IR)=0. 

RL(IR)=RL(IR) +RDLI (IR)#0T 

RDLR=ZDLF( IR)«( RL ( IR ) 4=AKRL ( IR ) -RLF ( IR) ) 

RLF( IR) =RLF ( I R) +RDLR 
IF(NORES.EQ.l )RLF ( IR)=0. 

RDH I ( IR )=R0HI ( IR ) +RDH( IR) ^OT 

IF( IMOOE.GE.O ) GO TO 245 
ROHK IR )=RDH( IR) 

RHFI(IR)=0. 

RFF I n I R)=( RDH I ( IR )^GNRS( IR ) )/0MGRE2( IR ) 

245 CONTINUE 

RHF( IR )=RDHI( IR ) =!^GNRS ( IR ) -TOMGRE ( IR) -RHF I ( IR) 

=!' -0MGRE2I IR)^RHFI I( IR) 

RHFI ( IR )=RHFI ( IR ) +RHF( IR) *DT 
RFFIK IR) = RhFII( IR)+RHFI( IR )^OT 

ROTATIONAL POSITION LIMIT 

RSF ( IR) =RLF( IR) +RHF( IR) +RIC ( IR) 

RSFX( IR)=RSF( IR) 

I F(ABS ( RSFX ( IR ) ) .GT.RPCSLI IR) ) RSFXI IR ) =S I G M RPOS L ( IR) 
*,RSFX( IR) ) 

ROTATIONAL POSITION COMMAND OUTPUT 
RSD( IR)=RSF(IR)+ROH( I R) /OMG RQ ( I R ) 

I F(ABS (RSD( IR ) ) .GT.RPOSLI IR ) ) RS D I I R ) =S I GN ( RPOSLI IR) 
-,RSD( IR) ) 

250 CONTINUE 

CPHLF=XCOS(RLF(l) ) 

SPHLF=XS INI RLF( 1 ) ) 

CTHLF=XC0S(RLF(2) ) 

STHLF=XSIN( RLF ( 2 ) ) 

C?SLF=XC0S( RLFI 3) ) 

SPSLF=XSIN( RLFI3) ) 

CPHSF = XCOS( RSFI 1 ) ) 

SPHSF = XSIN( RSFI I ) ) 

CTHSF = XCOSIRSFI 2) ) 

SThSF=XSINI RSFI2 ) ) 

CPSSF = XC0SIRSFI3) ) 

SPSSF=XSINIRSFI3) ) 

TRANSLATIONAL (CAB TO LINEAR DRIVE) 

transformation 

MATRI X 



TMTI 1,1 ) =CPSSF^CTHSF 
TMTI 2, 1 )=SPSSF 
TMTI3,1 )=-CPSSF*STHSF 

TMTI 1,2 )=-CPHSF*SPSSF*CTHSF + SPHSF*STHSF 
TMTI2, 2 )=CPHSF=i'CPSSF 

TMTI 3,2 )=CPHSF*SPSSF*STHSF + SPHSF=!‘CTHSF 
TMTI 1 ,3 )=SPHSF*SPSSF#CTHSF+CPHSF=icSTHSF 
T.MTI 2, 3 ) = -SPHSF*CPSSF 

TMTI 3, 3 )=-S?HSF*SPSSF*STHSF + CPHSF=f'CTHSF 
FNCLI 1 ) =STHLF#CPSLF 

FNCLI 2)=-CPHLF=«^STHLF=l^SPSLF-SPHLF*CTHLF 
FNCLI 3 )=SPHLF=!'STHLF=^'SPSLF-CPHLF=i'CTHLF 

PCHII 1) =FKNI 1 )-FNCLIl ) /AKXL 
FCHII 2) =FKNI2 ) -FNCL I 2) /AXYL 
FCHII 3) =FKNI3) 
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c 

c 



c 

c 

c 

c 

c 



c 

c 

c 



DO 270 IT = l ,3 

IF ( IMOOE.GE .0 ) GO TO 267 
FNCHK I T)=0. 

FNCHIK IT)=(FCHI< IT)*GNH( IT) ) /CMGT2( IT) 

267 CONTINUE 

FNCH( IT )=GNH( I T ) *FCH I ( I T) -TOMZT ( IT)*FNCHI( IT)- 

* CMGT2( IT I^S'FNCHI I ( IT) 

FNCHK IT)=FNCHI( IT)+FNCH( IT)*0T 
FNCHIK IT) = FNCHI K IT)+FNCHK IT)*0T 
T00H( IT )=G*(FNCH< IT)+FNCL( IT ) ) 

270 CONTINUE 

DO 300 IT=l,3 

TD0A( IT)=0. 

DO 280 JT=l,3 

TCDA( IT )=TDDA( IT ) +TDD H( JT) *TMT( IT,JT) 

230 CONTINUE 

TDDT( IT)=(TDDA( IT)-TDDG(IT) )*FN0T(IT) 

, DUMMY FOLLOW-UPS 

IF( IDUMFU.EQ.l) FOLUPS(IT)=TSD(IT) 

CHECK OPERATING AREA PARABOLIC LIMIT? 

IFIFOLUPSI IT) =i=TSI L( IT) ) 233,2 83,28 2 
292 ZTEMP=TSIL( IT)«TSIL( IT)*TMS (IT) 

IF( ( TPOSL ( I T)-ABS (TSD ( IT) ) ) .GE. ZTEMP) GO TO 283 
KFLAGI I T ) = l 

TDDT( IT)=-SIGN(ADDL(IT) ,F0LUPS ( IT) ) 

GO TO 294 

283 I F(KFLAG( IT ) .EQ.O ) GO TO 294 

IF(F0LUPS( IT)*TDDT( IT)) 236,286,285 

285 TDDT( IT )=0. 

GO TO 294 

286 KFLAG( I T)=0 

294 CONTINUE 

IF(A8S(TD0T(IT) ).GT.ADDL( IT)) TDD T ( I T ) =S I GN ( ADOL ( IT) 
*,TD0T( IT) ) 

IF( IMODE) 295, 297, 297 

295 CONTINUE 
TSIL( IT)=0. 

TDCSP( IT)=0. 

TCCSI ( I T)=0. 

TDCSI I ( IT)=TDDT( IT)/0MGTE4( IT) 

297 CONTINUE 

TDDS( IT)=TDDT( IT)-FOMGTE( IT)*TSIL( IT)-S0MGT2I IT)* 

* TDCSP( IT) -F0MGT3( IT)*TDC5I ( IT) -0MGTE4( IT)*TOCSI I ( IT ) 
TSIL( IT)=TSIL( IT)+TODS( IT)*OT 

TDCSP( IT) =TDCSP( I T ) +TS I L ( I T ) *0T 
TDCSK IT)=TDCSK IT )+TDCSP ( I T) *0T 
TDCSI KIT )=TDCS I I ( IT) + TDCS I ( IT)*DT 

TRANSLATIONAL POSITION COMMAND OUTPUT 
AND LEAD COMPENSATION 

TSD( IT )=TDOS( IT)*OMGQSQ(IT)+TSIL( IT)*GLEAD( IT)+ 

* TDCSPl IT ) + TI C( IT) 

3C0 CONTINUE 

60 CONTINUE 

RETURN 

END 



81 



ono 



c 

c 

c 



c 



c 

c 



c 

c 



initialize the filter parameters for 
THE drive program (SUBROUTINE SOIMOT) 



SUBROUTINE SOIMIC(DTP) 

PROGRAMMED BY 0. ASTILL 



(CSC) 



12 SEPT 1975 



COMMON/ WSHOUT/AKR (3) , AKl ( 3 ) , OMGR ( 3 ) ,OMGT (3 ) , ZETR( 3 ) 
* t ZETT ( 3 ) 




R( 3) 



DO 200 IN=1,3 



AKRL( IN)=OMGRL( IN)/(2.=i'ZETRL(IN)) 

SCMGR( IN)=OMGR( IN 

ZCMGRd N) = 2 .’■-OMGR ( IN)«ZETR( IN) 

OMGQSQ( IN ) = U/OMGTQ( IN)**2 
GL£AD( IN)=2 .*ZETAQ( I N ) /OMGT Q ( I N ) 

TGMGRE ( IM) =2.-0MGRE( IN ) 

0MGR£2( IN)=OMGRE( IN)-*2 

0MGT2( IN)=OMGT( IN)^*2 

TCMZT( I N) = 2 .*OMGT ( IN)=!'ZETT( IN) 

FCHGTE( IN )=A.=«'OMGTE( IN ) 

SCMGT2( IN)=6.*0MGTE( IN)**2 
FCMGT3( IN)=4.*0MGTE( IN)**3 
CMGTE4( IN)=OMGTE( IN)*#4 

ZDLF( IN> = 1.-EXP(-0T?*2.*ZETRL( IN)*QMGRl ( IN) ) 



2 CO CONTINUE 



RETURN 

END 



C BLOCK DATA SUBPROGRAM TO INITIALIZE THE NON-WASHOUT 
C FILTER PARAMETERS 
C 

BLOCK DATA 
r 

CCMMON/MOTION/OMGRL( 3) ,OMGPQ( 3) ,0MGRE(3) ,0MGTE(3) 

I , ZETAQ(3 ) , ZETRLI3 ) ,ZOLF( 3 ) ,GNRS (3) ,AKRL ( 3) ,SGMGR (3) 

2, ZCMGR(3),GMGQSQ(3) ,T0MGREI3) ,OMGRE2(3) ,CMGT2(3) 

3, TOMZK 3) tFOMGTE ( 3) , S0MGT2( 3) , FCMGT3(3) ,0MGT£4(3) 

4, 0MGTG(3 ) ,GLEAO( 3) 

C 

DATA ZETRL/3*.707/t ZETAQ/.7,.?, 50./ 

DATA GNRS/3=i'l ./,0MGRL/3=i'l .2/,OMGRQ/3.5, 7.,8.5/ 

DATA OMGTQ/4. 1,7. , I 30. / ,OMGTE /3=i=. 0 5/ 

DATA 0MGRE/3=«= .01/ 

C 

END 
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c 



GRMS8Y HUMAN DYNAMIC ORIENTATION MODEL SUBROUTINE 



C 

C 




C 

C 



10 

15 



C 

C 



20 

25 



)tWbrU(JJjrNl3j fYl3J 
DIMENSION W(N,3) tSF(N,3) 

DATA DT/1 .0/,NITP/10/ 

REWIND lOUT 
REWIND 4 

NOT=NSP=f'DT/NlTP 
INDT=(NDT-2)/2 
WRITE! I OUT, AGO )0T, I NOT 

CANAL SPECIFICATIONS. 



DO 10 1=1, A 
READ(A, 530) 
CCNT INUE 
DO 15 1=1,3 
READ! A, 530) 
CONTINUE 
READ(A, 530) 
READ! A, 530) 
READ! A, 550) 



TCI I , 1 ) ,TC! I ,2 ) , TC ! I, 3) ,TC! I f A ) 



TPC { I , 1 ) ,TPC( I , 2) ,TPC (I , 3) , OVC ( I ) 

CC( 1 ) ,CC! 2) ,CC! 3) ,CC( A) 

GKC ( 1 ) , GKC ( 2 ) , GKC ! 3 ) , GKC ( A ) 
FSCC,TSCC,SSCC 



27 

C 

C 



32 



36 



A5 

50 



CALL EULER(FSCC,TSCC,SSCC,CTC) 

OTOLITH SPECIFICATIONS. 

DO 20 1=1,3 

READ! A, 550) TO ! I , 1 ) , TO ! I , 2 ) , TO ! I , 3 ) 

CCNT INUE 
DO 25 1=1,2 

READ! A, 550) TPO! I , 1 ) , TPO! I , 2 ) , D VC ! I ) 

CONTINUE 

READ! A, 550) cod ) ,C0!2 ) ,C0!3 ) 

READ! A, 570) GKO! 1) ,GK0(2) , GKO ! 3) ,GKS( 1) ,GKS!2) ,GKS!3) 
READ! A, 570) FOTO , TOTO , SOTO , S AC F AC, OSPG , DFAC 
CALL EULER! COTO, TOTO, SOTO, CTO) 

DO 27 1=1,3 
CS!I)=CO( I )-SACFAC 
CONTINUE 

INITIALIZATION 
DO 32 1=1,4 

READ! A, 556)XCH( I ) ,YCH! I ) , ZCH! I ) 

CONTINUE 

A!1)=XCH(A) 

A! 2)=YCH! A ) 

A! 3)=ZCH! A) 

CALL COTRN! A,CTC, 1 ,W0) 

DC 36 1=1,3 

READ!A, 550) XOH! I ) , YOH ! I ) , ZOH ! I ) 

CONTINUE 
A! 1)=X0H!3) 

A! 2)=YQH! 3 ) 

A! 3 )=Z0H<3 ) 

CALL COTRN! A, CTO, 1 ,A0) 

DO A5 1=1,3 

R6AO(A,5AO)XC! I) ,YC(I ) ,ZCd ) ,WNCO( I) ,WMCL( I ) 

CONTINUE 
DC 50 1=1,2 

READ(A, 550) XO (I ) , YO (I ) , ZO ! I ) 

CONT INUE 



REAC(A, 550 ) DOLDd ),D0LD(2 ) ,D0 lD(3) 
READ! A, 550) WSFOI I ) , WS FO ! 2 ) , WS FO ! 3 ) 
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55 



REA0(4,550)W0F(1 ) , WOF ( 2 ) , WOF ( 3 ) 

R EAO( 4, 580 ) TOVFL , TOP OS ,TNC »FN0I SE 
FP=1 .0-EXP( -DT/TOPOS) 

FD=T0VEL 

FN( 1) = EXP(-0T/TNC) 

FN( 2)=TNC=!'( l.-FN( 1 ) )/0T-FN( 1) 

FN( 3) = l.-TNC=<'( 1 .-PN( 1 ) ) /DT 

MAIN PROGRAM CYCLE. 



COMPUTE CURRENT STIMULUS IN HEAD COORDINATES 
(EVERY DT/NITP SEC.) 

1. ANGULAR ROTATION VECTOR (TWH) AT (TIME). 

2. SPECIFIC FORCE VECTOR ( TAH ) AT (TIME+DT/2) 

3. TRUE DOWN VECTOR* AT (TIME+OT/2). 

00 450 ITIME=1 , I NOT 

DO 100 I=1,NITP 

T IME=( I TIME -1 >*0T + I*DT/N ITP 



J = 20.=!'TIME+1.5 
K=J+10 



TAH( 1)=SF( K,l)/32.2 
TAH(2)=SF(K,2)/32 .2 
TAH(3)=SF(K,3)/32.2 
TVsH( 1 ) = W( J, 1) 

TV.H(2 ) = W( J, 2) 

TV/H( 3) = W( J, 3) 

transform to sensor COORDINATES. 

CALL COTRN( TWH,CTC,0,TWS) 

CALL COTRN( TAH,CTQ,0t TAS ) 

SENSOR STIMULATION (EVERY DT/NITP SEC.): 

USING CURRENT STIMULUS VALUES, UPDATE STATE VECTORS 
FOR 3 CANALS ( XC , YC, AND ZC) , AND 3 OTOLITHS 
(XO, YO, ZO), AND COMPUTE AFFERENT FIRING RATES 
(CSX, CSY, CSZ). 

S = TWS ( 1 ) 

CALL SVUPO(XC,TPC,DVC,S,CSX ,CC, 3,4) 

S = T'/^S( 2 ) 

CALL SVUPD( YC,TPC,DVC,S,CSY ,CC, 3,4) 

S=TWS(3 ) 

CALL SVUPD( ZC,TPC,DVC,S,CSZ,CC,3,4) 

S=TAS( 1 ) 

CALL SVUPO( XO,TPO,DVO,S,OSX,CO, 2,3) 

S=TAS(2) 

CALL SVUPD( YO,TPO,OVC,S,OSY,CO, 2,3) 

S = TAS( 3 ) 

100 CALL SVUP0( ZO,TPO,DVQ,S,OSZ,CS, 2,3) 

optimal ESTIMATOR (UPDATE EVERY DT SEC.): 

GET CANAL AND OTOLITH SYSTEM, STATE ESTIMATES FROM 
STEADY STATE KALMAN FILTERS. 

CALL SSKF(XCH,CSX,TC,CC,GKC,4) 

CALL SSKF( YCH,CSY ,TC,CC ,GKC ,4) 

CALL SSKF( ZCH,CSZ ,TC,CC,GKC,4) 

CALL SSKF(X0H,0SX,T0,C0,GK0,3) 

CALL SSKF( YOH,OSY ,T0,CC,GK0,3) 

CALL SSKF( Z0H,0SZ ,TO,CS,GKS ,3) 

ENTER ROTATION RATE ESTIMATE VECTOR (CANAL ESTIMATE). 
EVS( 1 ) =XCH(4) 

EUS( 2)= YCH( 4) 

E'aS(3 ) -ZCH( 4) 

ENTER SPECIFIC FORCE ESTIMATE VECTOR (OTOLITH ESTIMATE) 
EAS( 1 )=XOH( 3) 

EAS( 2) = Y0H( 3) 
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c 

c 

c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 

c 

c 



c 

c 

c 



SACCULE NON-LINEARITY. - 
E AS (3 ) =A(MAX1 ( .6* ( ZOH( 3 ) +. A1 69 ) - . 4 1 69 1 - . 4 1 69 ) 

otolith ESTIMATE TO VALUE HELD 
®|£9^ECONSI DERATION OF SACCULE NON-LINEARITY. 

ThE^ofoLITH StTmATe")^ EFFECTS ONLY DIRECTION OF 

CALL NORM(EAStY) 

DC 130 1=1,3 

130 EAS(I) = Y(IJ *SQRT ( XOH ( 3) *=!=2 + Y0H( 3 )**2 + Z0H (3 ) *=*'2) 

TRANSFORM TO HEAD COORDINATES 
CALL COTRN( EWS,CTC,l ,EWH) 

CALL COTRNI EAStCTO, 1,EAH) 

DOWN AND W ESTIMATOR lUPDATE EVERY DT SECONDS). 
COMBINE OTOLITH AND CANAL ESTIMATES TO FORM NEW 
ESTIMATE OF: 

1. PERCEIVED DOWN (DNEW) AT (TIME+OT/2). 

2. PERCEIVED ACCELERATION (ACC) AT (TIME+DT/2). 

3. PERCEIVED ANGULAR VELOCITY (WTOT) AT (TIME). 

^40 CALL DCWN(DOLD,EWH,EAH,AO,WSFO,FD,CT,TOPOS,DFAC,WOF 
WNCO ,WNCL ,FN, lOUT) 

450 CONTINUE 

FORMAT STATEMENTS. 

400 FORMAT! F6. 3, 14) 

52C FORMAT( E15.8, 15 ) 

5 50 FORMAT! 4El 5.8) 

540 FORMAT! 5E15.8) 

550 FORMAT! 3E1 5.8) 

570 FORMAT! 6E1 2.5) 

580 FORMAT! 4E15 .3 ) 

RETURN 

END 



C DOWN ESTIMATOR AND W ESTIMATOR 
C 

SUBROUTINE DOWN ! OOLD , W N ,SN, SO, WSFO ,T , DT , TOP S, OFAC, WOF 
WNCO ,WNL, FN, lOUT ) 

C 

C DOWN IS DETERMINED BY RELYING ON LOW FREQUENCY OTOLITH 
C ESTIMATES, CANAL ESTIMATES THAT ARE CONSISTENT WITH 
C HIGH FREQUENCY OTOLITH ESTIMATES, AND THE HIGH FREQUENCY 
C PORTION OF CANAL ESTIMATES NOT CONFIRMED BY OTOLITHS. 

C 

C W IS DETERMINED BY CANAL ESTIMATES PARELLEL TO DOWN, 

C ROTATION RATE OF DOWN, AND HIGH FREQUENCY PORTION OF 
C CANAL ESTIMATES PERPENDICULAR TO DOWN MINUS ROTATION 
C RATE OF DOWN. 

C 

DIMENSION DOLD! 3 ) ,WN ! 3 ) , SN( 3 ) , SO ( 3 ) , WSFO ! 3 ) , F ( 3 ) 

1 , X !3 ) , WOF! 3 ) ,RSCC ( 3) , ROTO ! 3 ) , RTOT ! 3 ) , DNE W ( 3 ) ,R?OS ( 3) 

2, WPERP ( 3 ) ,D AVG!3) ,WPARE (3) , WTOT ! 3) , ACC ( 3 ) , ANG! 3 ) 

3, WNCO! 3 ) , WNL! 3 > , FN! 3 ) ,WNC!3) , WNCH(3) , HRO TO ! 3 ) , WSF ( 3 ) 

4, WCD(3) ,W0DN(3) 

C 

SFMAG=SQRT!SN! 1 ) =<^SN ! 1 ) +SN ( 2 ) *SN ! 2 > +SN ( 3 ) *SN ( 3 ) ) 
FPCS=1. O-EXP!-! (SFMAG/ OFAC )=<'-( .2 5 ) )*DT/TOPS ) 

TCVEL=T 

F ! 1 ) = EXP(-OT/TDVEL ) 
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F( 2)=TDVEL*(1.-F( 1) )/DT-F( 1) 

F(3)=1.-T0VEL=»( 1 .-F(l ) )/0T 
CALL CROSS ( SO, SN, WSF) 

CALL NORM(WSF,X) 

CALL VANG( SO,SN, ANGSF) 

OG 10 1=1,3 
WSF(I)=ANGSF*X( I) 

WGF( I) =F( 1 )*WOF( I )+F(2 )*WSFO{ I ) +F(3)=!»WSF( I ) 

WCO( I )=WSF( I)-WOF( I ) 

10 WSFO( I )=WSF ( I ) 

WOOM=SQRT( WGO ( I )=*WOD( 1 ) +W00( 2 ) *W00 ( 2 ) +W00( 3 )*WOn( 3) ) 
CALL NORM! WOO, WOON) 

WCPWO=WN ( 1 ) =SWOON ( 1 )+WN ( 2 ) *WOON ( 2 ) +WN ( 3 ) *WGON ( 3 ) 
IF(WCPWD) 12,12,11 

11 WCPWD=0.0 

12 WMAG=-WCPWO*OT 
00 13 1=1,3 

X( I)=W00N( I )*AMIN1(WMAG,W00M) 

WNC( I )=-WN( I )*0T-X( I ) 

WNL( I) = FN( 1 )=^WNL( I ) + FN( 2)=i'WNC0( I )+FN( 3)*WNC(I ) 

WNCO( I) =WNC( I ) 

WNCH( I ) =WNC ( I )-WNL ( I ) 

13 X ( I ) = X( I )+WNCH( I ) 

CALL CROSS (OOLO, SN, ROTO) 

CALL NCRM( RQTO,F) 

WGPARM=V')OF { D^^'Fd ) +WOF ( 2 ) *F ( 2 ) +WOF ( 3 ) ( 3) 

OC 15 I =1 , 3 

ROTO( I )=WOPARM*F ( I ) 

15 HRCTOl I )=ROTO( I ) /2. 

CALL ROTATE (OOLO ,HROTO ,F) 
WCPARM=X(1)*F(L)+X(2)*F(2)+X(3)*F<3) 

DC 20 1=1,3 

RSCC( I) = (X( I)-WCPARM*F( I) ) 

20 RTOT( I )=RSCC( I)+ROTQ( n 

CALL RCTATE1D0LD,RT0T,DNEW) 

CALL VANG( ONEW, SN ,FEE ) 

OUT — POnc^PPc: 

CALL CROSSIDNEW, SN,RP0S) 

CALL NORM(RPOS,X) 

OC 30 1=1,3 
30 RPOS( I)=PHI*X( I ) 

CALL RQTATE(DNEW,RPOS,X) 

CALL NORM(X,ONEW) 

DO 40 1=1,3 

40 X( I )=DOLO( I )+ONEW ( I ) 

CALL NORM(X,DAVG) 

WPARM=WN( 1 ) ^OAVG( 1 ) +WN ( 2 ) *0 AVG ( 2 ) +WN ( 3 ) *OAV G( 3) 

CALL CROSS(ONEW,OOLD,WPERP) 

CALL NORM! WPERP ,X) 

CALL VANGIOCLD, ONEW, PHI ) 

DO 50 1=1,3 

WPERP ( I ) = X( D^i'PHI/DT 

WPARE ( I ) =WPARM=!=OAVG( I ) 

DCLD( I ) =ONEW( I ) 

SG( I )=SN( I ) 

WTQT( I ) =WPERP ( I )+ WPARE ( I ) 

ACC( I ) = OFAC*ONEW( I )-SN( I ) 

50 CONTINUE 

45 WRITE! I OUT, 100) ACC ( 1) , ACC ( 2 ) , ACC ( 3 ) , WTOT ( 1 ) , WTOT ( 2 ) 
=!=, WTOTI 3 ) ,ONEW( I ) ,ONEW( 2) ,DNEW( 3) 

ICC FCPNAT! 9E14.7) 

RETURN 

END 
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c 

c 

c 

C 

c 

c 



c 



c 



c 



c 



^9 0E\/EL0P a cost algorithm based on the 

AIRCRAFT^MOTION^^*^ PERCEPTION OF SIMULATOR AND 



SUBROUTINE COST ( C AO , CAVGt CRMS ) 

COMMON/ FL I GHT/XAI 4884 t3 ) t XR ( 4884 » 3 ) , P , C , I P AR 

01 MENS ION AA(3) tAS(3) ,WA( 3) ,WS( 3),DWNA( 3 ),0WNS( 3) 
DATA WTH,ATH/. 0146, .005/ 

CAVG=0. 

CRMS=0. 

CAC=0. 



CV<AV=0. 

CAAV=0. 

CCAV=0. 



C^RMS=0 . 
CARMS=0 . 
CCRMS=0. 



A0V^=0. 

ACA=0. 

ACC=0. 



READ IN THE DATA FROM THE TWO FILES 

REWIND 9 
REWIND 8 



REAC(9, 30)DTS,NSS 
REAO( 8, 30 ) DTA,NSA 

COMPUTATION OF THE COST ESTIMATION VALUES 



DO 40 1=1, NSA 
C 

R£AD( 9, 101 AS( 1 ) , AS ( 2 ) , AS ( 3 ) , WS ( 1 ) , WS ( 2 ) , WS ( 3 ) ,DWNS( 1 ) 
*,DWNS( 2 ) ,DWNS( 3) 

REAC(3, 10)AA(1 ),AA(2) ,AA(3) ,WA( 1 ) ,WA(2) ,WA( 3) ,0WNA(1) 
=!=,DWNA( 2 ) ,DWNA( 3) 

C 

CW1=WA( 1 )-WS( 1 ) 

0W2 = WA( 2 )-WS( 2 ) 

0W3=WA( 3)-WS( 3) 

C 

D A1=AA( 1)-AS( 1 ) 

CA2=AA( 2)-AS(2) 

CA3=AA( 3)-AS(3) 

C 

CD1=DWNA( 1 ) -OWNS ( I ) 

CC2=0WNA( 2)-DWNS{ 2) 

CD3=DWNA( 3) -DWNS( 3) 

C 

DW=SQRT {0WI*DW1+DW2-DW2+DW3*DW3 ) 

CA=SQRT (OAl*OAl+DA2-OA2+DA3*OA3 ) 

DD = SCRT (DO 1*001 +DD2’:'0D2 + D03=!--003) 

C 

AWA=SQRT( WA { 1 ) *WA ( 1 ) +WA( 2 )*WA( 2 ) +W A ( 3 ) *W A ( 3 ) ) 
AAA=SQRT( AA(n*AA(l)+AA(Z)*AA(2)+AA(3)*AA(3)) 

C 

C'/v = OW/AMAXl ( AWA, WTH) 

CA=CA/AMAX1 ( AAA, ATH) 

CC = CO 

r 

CWAV=CWAV+CW 
CAAV=CA AV+CA 
CCAV=CDAV+CD 
C 

CWRMS=CWRMS+CW*CW 
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C;iRMS=C ARMS+CA*CA 
CCRMS=CORMS+CO=i'CD 
C 

ACU=ADW+OW 

ACA=ADA+OA 

ACC=A00+00 

C 

40 CONTINUE 
C 

CWAV=CWAV/NSA 

CAAV=CAAV/NSA 

CCAV=COAV/NSA 

CAVG=CWAV+CAAV+COAV 

C 

CWRMS=SQRT (CWRMS/NSA ) 

CARMS=SQRT( CARMS/NSA) 
CDRMS=SQRT(CDRMS/NSA) 
CR.«S=CWRMS+CARMS+CORMS 
C 

ACw=AOVj/NSA 
ACA=ACA/NSA 
AOO=AOD/NSA 
C AD=AOW+AOA+AOD 
C 

V«RITE(6,70) AOW, ADA, ADD, CAD 
70 FCRMAK//’ AOW = • , F 1 1 . 6 , ’ A0A=',F11.6, 

FI 1 .6, • TQTAL=' ,F10.6) 

^.RITE(6 ,50)CWAV,CAAV,C0AV,CAVG 
50 FCRMATt/’ CWAV = • , FIO .6 , • CAAV=',F10.6 

FIO. 6,' TOTAL=' ,F10.6 ) 

WRITE(6 ,60)CWRMS, CARMStCDRMStCRMS 
60 FCRMAK/' CWRMS = ’ , F9 .6 , • CARMS=’,F9.6 

=f=,F9.6, • TCTAL=' ,910. 6/) 

30 FCRMAK F6.3 , 14 ) 

10 FCRMAK9E14.7) 

C 

RETURN 

END 



ADD = « 

• CDAV= 

• CDRMS 
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c 

c 

c 

c 



50 

20 

30 

60 

40 

10 

C 



SU6RCLTINE TO DIFFERENTIATE THE TIME SERIES 
CONSISTING OF THE ROTATIONAL POSITION 
COMMANDS FROM THE DRIVE PROGRAM 

SLBRCUTINE 01 FF ( DT , Y, YD , NSP ,N ) 

DIMENSION Y(N,3) »Y0(N,3) 

DC 10 1=1, NSP 

DO 40 J=l,3 

I FU-I ) 20,20, 50 

IF( I-NSP) 60,30,30 

YC( I,J)=(Y{ I + l, J )-Y{ I,J))/OT 

GO TO 40 

YC(I,J) = (Y(I,J)-YU-1,J))/0T 
GC TO 40 

YC( I, J)=(Y( I+l , J)-Y( I-1,J) ) /(2.*DT) 

CONTINUE 

CENT INUE 

RETURN 

ENC 



C SLBRCUTINE TO TRANSFORM THE COORDINATE SYSTEM OF A VECTOR 
C 

SUEROUTINE COTRN ( A , B , N ,C ) 

C 

C CCCRDINATE TRANSFORM: 

C 

C FRCM HEAD TO SENSOR IF N=0 (DIRECT TRANSFORMATION) 

C FRCM SENSOR TO HEAD IF M=-l (INVERSE TRANSFORMATION) 

C FRCM SENSOR TO HEAD IF N=1 (INVERSE TRANS FORMAT ION 

C USING THE TRANSPOSE FOR ORTHOGONAL MATRICES) 

C 

C A=CRIGINAL VECTOR 

C B=TRANSFORMATION MATRIX 

C C=TRANSFORMED VECTOR 

C 

DIMENSION A{3),B{3,3),C(3),BI(3,3),L(3),M(3) 

C 

IF(M20,10,30 

C 

10 DC 15 1=1,3 

C(I)=B( I,1)*A( 1) + B(I,2)^A(2) + B( I ,3)-A(3) 

15 CONTINUE 

RETURN 
C 

20 CALL MCPY(3,6I , 3, 3,0) 

CALL MI NV( B I ,3 ,D, L,M) 

DO 35 1=1,3 

C(I)=BI(I,1)-A(1)+8I(I,2)=!=A(2)+BI(I,3)*A(3) 

35 CONTINUE 

RETURN 
C 

30 DC 40 1=1,3 

C { I )=B( 1 , I )*A( 1 )+8( 2, I ) -A( 2 )+6( 3, 1 )*A( 3 ) 

40 CONTINUE 

C 

RETURN 

END 
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SIBROUTINE EULER( F,T,S,CT) 

PRGCUCE DIRECTION COSINE MATRIX (CT) GIVEN EULER ANGLES 
(F,T, AND S) 

DIMENSION CT{3,3) 

CT( It 1 )=COS(S)=!'COS(F )-COS (T)*SIN(F)^SIN(S) 

CT( 2, 1)=-(S IN(S)*COS(F)+COS(T)*SIN{F)*CGS(S)) 

CT(3,1 )=SIN(T)=!'SIN(F) 

CT(lt2)=C0S(S)*SIN(F )+COS ( T ) *COS ( F )*S IN ( S) 

CT ( 2 » 2 ) =C OS (T ) *C0 S { F ) =^^0 OS ( S ) - S I N ( S ) *S I N ( F ) 
CT(3,2)=-SIN(T)*C0S(F) 

CT{ 1,3) = SIN(T)-SIN(S) 

CT(2,3 )=SIN(T)*C0S(S) 

CTO, 3) =C0S(T) 

RETURN 

END 



SLBROUTINE CROSS (A, 8,0 
C = A X B 

DIMENSION A(3) ,B( 3) ,C( 3> 
C(i)=A(2)-6(3)-A(3)*B(2) 
C ( 2) = A{ 3)=!'B(1 )-A( 1)*3(3) 
C(3)=A( 1)«B(2)-A( 2)-B( L > 

RETURN 

ENC 



SLBROUTINE V ANG ( A , B , PHI ) 

°HI = ANGLE BETWEEN A AND 8. 

DIMENSION A(3),B(3),AN(3),3N(3) 

CALL NORM! A, AN) 

PhI = ARCOS ( ANO ) '-^BN { 1 ) +AN< 2 ) *6N { 2 ) +AN ( 3 ) - BN< 3) ) 

RETURN 

ENC 



C 



SUBROUTINE NORM! A, AN) 

AN = UNIT VECTOR IN DIRECTION OF VECTOR A. 

DIMENSION A(3),ANOI ' 

AM = SCRT( A( I )*A( 1 ) + A( 2 I'i'Al 2) + A( 5 ) *A(5) ) 

I F ( AM~ 1 . E~0 6 ) 5,5,6 

5 AM=l.E-06 

6 DC 10 1=1,3 
1C AMI)=A( I)/AM 



RETURN 

END 
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SLBROUT INE SVUPO < X , T , D, S t Y ,C , N » M ) 

ST/iTE VECTOR UPDATE: 

X(NEV.) =T#X(OLO) + O^S 
Y<NEW) = C=f=X(NEW) 

WHERE: N IS DIMENSION OF STATE VECTOR 
X IS STATE VECTOR 
T IS TRANSITION MATRIX 
D IS DRIVING VECTOR 
S IS STIMULUS 

Y IS OUTPUT (AFFERENT FIRING RATE) 

DIMENSION X(N),T(N,N)tO(N),C(M),R(9) 

DC 5 1=1, N 
= R( I ) = X( I) 

DC 10 1=1, N 
X(I) = D( I)=!‘S 
DC 10 J=1,N 

10 X(I)=X( I)+T(I,J)«R(J) 

Y=C(M)*S 
DC 20 I=1,N 
20 Y = V+C( I )*X( I) 

RETURN 

END 



SUBROUTINE SSKF ( XH , Y , TM ,C , GK, N ) 

STEACV STATE KALMAN FILTER (UPDATE EVERY DT SECONDS) 

WHERE: XH IS STATE VECTOR ESTIMATE 
TM IS TRANSITION MATRIX 
GK IS KALMAN GAIN MATRIX 
Y IS SENSOR SYSTEM OUTPUT 
C IS OUTPUT MATRIX 

DIMENSION XH(N),TM(N,N),C(N),GK(N),S(9) 

DO 40 I=1,N 
S( I)=0.0 
DC 40 J=1,N 

40 S ( I ) = S( I )+TM( I , J ) *XH( J) 

EM=0.0 
DC 45 1=1, N 
45 EM = EM+S ( I )«C( I ) 

DC 50 J=1,N 

50 XH(J)=S( J)+GK(J)^(Y-EM) 

RETURN 

END 



SLEROUTINE ROTATE (A,R,AR) 

AR = A ROTATED ABOUT R BY ANGLE (RAD) EQUAL TO THE 
MAGNITUDE OF R 

D I MENS ION A (3 ) ,R(3 ) , AR( 3) ,AP( 3) , APN( ?) 

CALL CROSS (R,A,AP) 

CALL NORM (AP,APN) 

AMAG=SQRT(A(1 j’SA(i)+A(2)*A( 2 ) +A ( 2 ) =^A ( 3 ) ) 
PHI=SQRT(R( 1)*R( I )+k(2)=«'R(2)+R{3)*R(3) ) 

00 10 1=1,3 

10 AR ( I )=AMAG-SIN( PH I )=!'APN( I )+COS (PHI )*A( I ) 

RETURN 

END 
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